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ABSTRACT 


This  thesis  describes  the  development  of  computer  algorithms  for  experimental 
measurement  control  and  subsequent  signal  processing  of  transient  signatures  to  syn¬ 
thesize  scattering  impulse  responses  of  scale  model  targets.  The  theories  behind  tran¬ 
sient  scattering  are  considered  in  order  to  construct  the  algorithms.  Up-to-date 
hardware  and  softv/are  technology  are  selectively  implemented  to  optimize  the  resultant 
signal-to-noise  ratio.  The  detailed  explanation  of  the  hardware  and  software  operation 
is  provided.  A  noise  model  for  the  system  is  also  discussed.  Measurement  and  pro¬ 
gramming  validations  are  described,  where  comparisons  are  made  with  numerical  com¬ 
putations  for  transient  scattering  by  selected  canonical  targets.  Results  are  documented 
for  futher  research  in  the  transient  scattering  problem. 


iii 


Accession  For 

NT IS  GRAAI 

art 

DTIC  TAB 

o 

Unannounced 

Just  If leal  Ion — 

□ 

By 

i 

Distribution/ 
Availability  Codes 
Avail  and/or 
Special 


1  Dlst 


TABLE  OF  CONTENTS 


I.  INTRODUCTION  . 1 

A.  BACKGROUND  . 1 

B.  OVERVIEW  OF  THIS  EFFORT . 2 

II.  SYSTEM  FACILITIES  AND  HARDWARE  . 4 

A.  LABORATORY  DESCRIPTION  . 4 

B.  ANECHOIC  CHAMBER  . 6 

C.  HP54120T  DIGITAL  PROCESSING  OSCILLOSCOPE  . 6 

D.  GURU  II  HARDWARE  . 8 

III.  THEORY  OF  TRANSIENT  SCATTERING  MEASUREMENTS  . 11 

A.  FREE-FIELD  RANGE  SYSTEM  REPRESENTATION  . 11 

B.  MATHEMATICAL  MODEL  . 13 

C.  TRANSIENT  RESPONSE  SOLUTION  . 14 

IV.  ACQUISITION  AND  SIGNAL  PROCESSING  ALGORITHMS  . 17 

A.  GURU  II  INTERFACE  . 17 

1.  GPIB  Principal  . 17 

2.  Basic  Software  Requirements  . 18 

B.  ACQUISITION  ALGORITHM  . 18 

C.  DECONVOLUTION  ALGORITHM  . 20 

D.  EXAMPLE  OUTPUT  . 25 

V.  EXPERIMENTAL  RESULTS  . 27 

A.  CANONICAL  MODEL  SCATTERING  . 27 

B.  SCALE  MODEL  SCATTERING  . 32 

C.  SYSTEM  NOISE  ESTIMATION  . 33 

VI.  CONCLUSIONS  . 43 

A.  SUMMARY  . 43 

B.  FUTURE  CONSIDERATIONS  . 44 

iv 


APPENDIX  A.  THE  IMPULSE  RESPONSE  DECONVOLUTION  . 46 

APPENDIX  B.  THE  ACQUISITION  QUICK  BASIC  PROGRAM  . 67 

APPENDIX  C.  CHARACTERISTICS  OF  LOW  LOSS  CABLE  . 78 

LIST  OF  REFERENCES  . 79 

INITIAL  DISTRIBUTION  LIST  . 81 


LIST  OF  TABLES 


Table  1.  IIP54120T  CHANNEL  CHARACTERISTICS  . 7 

Table  2.  HP54I20T  TDR  AND  TLMEBASE  CHARACTERISTICS  . 7 

Table  3.  TYPICAL  NOISE  .  8 

Table  4.  DATA  OUTPUT  TRANSFER  RATE  AND  DATA  RECORD 

LENGTHS  . 8 

Table  5.  FREQUENCY  DOMAIN  FUNCTION'S  . 11 

Table  6.  MATHEMATICAL  NOTATION  . 14 

Table  7.  CALIBRATION  SPHERE  PARAMETERS  . 23 

Tabic  8.  90  DEGREE  THICK  WIRE  PARAMETERS  . 26 

Table  9.  12.2  CM  SPHERE  PARAMETERS  . 28 

Table  10.  45  DEGREE  THIN  WIRE  PARAMETERS  . 29 

Tabic  11.  90  DEGREE  THIN  WIRE  PARAMETERS  . 30 

Tabic  12.  8.1  CM  SPHERE  PARAMETERS  . 31 

Table  13.  15.5  CM  SPHERE  PARAMETERS . 33 

Table  14.  DIMENSIONS  OF  SCALE  AIRCRAFT  MODEL  l  . 34 

Tabic  15.  DIMENSIONS  OF  SCALE  AIRCRAFT  MODEL  2 . 35 

Tabic  16.  DIMENSIONS  OF  SCALE  AIRCRAFT  MODEL  3 . 36 

Table  17.  DIMENSIONS  OF*  SCALE  AIRCRAFT  MODEL  4 . 37 

Table  18.  STATISTICS  OF  NOISE  WAVEFORM  . 39 

Table  19.  STATISTICS  OF  NOISE  MODEL  . 40 

Table  20.  STA  TISTICS  OF  THICK  WIRE  SIGNAL  PLUS  NOISE  WAVEFORM  41 

Table  21.  CHARACTERISTICS  OF  PORT  RETURN  CABLES  . 78 

Table  22.  INSERTION  LOSS  . 78 


VI 


LIST  OF  FIGURES 


Figure  I.  General  layout  of  the  scattering  range  . 5 

Figure  2.  System  representation  of  the  TESL  . 12 

Figure  3.  Deconvolution  input  waveform  . 21 

Figure  4.  Overlay  . 23 

Figure  5.  Subtraction  waveform  . 24 

Figure  6.  Thick  wire  measurement  . 25 

Figure  7.  The  12.2  cm  diameter  sphere  validation . 28 

Figure  8.  The  10  cm  thin  wire  validation  at  45  degrees . 29 

Figure  9.  The  10  cm  thin  wire  validation  at  90  degrees . 30 

Figure  10.  The  8.1  cm  diameter  sphere  validation  . 31 

Figure  11.  The  15.5  cm  diameter  sphere  validation  . 32 

Figure  12.  Vertical  and  horizontal  wing  scale  aircraft  model  1  34 

Figure  13.  Scale  aircraft  model  1  at  000,  030,  090,  180  degree  aspect  angle  . 35 

Figure  14.  Scale  aircraft  model  comparisons  at  000  degree  aspect  angle  . 36 

Figure  15.  Scale  aircraft  model  comparisons  090  degree  aspect  angle  . 37 

Figure  16.  Noise  waveform  statistics  . 39 

Figure  17.  Noise  power  . 40 

Figure  18.  Signal  plus  noise  waveform  statistics  . 41 


vii 


ACKNOWLEDGEMENTS 


I  wish  to  thank  Professor  Michael  A.  Morgan  for  his  thoughtful  and  patient  guid¬ 
ance  throughout  the  course  of  this  research. 

This  thesis  is  dedicated  to  my  child  and  my  wife^M^who  for  ovei  two  and  a  half 
years  has  kept  the  home  fires  burning  while  I  have  pursued  my  education.  I  could  not 
have  done  it  without  her. 


I.  INTRODUCTION 


A.  BACKGROUND 

In  circuit  theory,  it  is  well  known  that  the  characteristics  of  a  linear  time-invariant 
system  can  be  obtained  by  exciting  the  circuit  with  an  impulse  function.  The  theory  of 
transient  scattering  measurements  is  analogous  to  this  method,  and  it  also  provides  a 
unique  parametric  description  of  the  scatterer  via  the  complex  pole  locations  in  the 
Laplace  domain  transfer  function.  There  are  two  approaches  to  obtaining  the  impulse 
response.  One  utilizes  frequency  domain  methods  and  the  other  employs  time  domain 
techniques.  The  former  has  to  sweep  frequencies  over  the  entire  band  but  the  latter  uses 
direct  pulse  excitation.  As  a  result,  the  latter  can  save  time  and  expense.  Therefore, 
direct  transient  scattering  measurements  are  a  viable  alternative  to  the  more  conven¬ 
tional  CW  method  for  obtaining  impulse  responses  of  targets. 

In  1965  Kennaugh  and  Moffatt  [Ref.  1J  investigated  the  properties  of  the  impulse 
response  waveforms  backscattered  from  metallic  targets.  The  approximations  for  tran¬ 
sient  scattering  and  impulse  responses  they  made  formed  the  basis  for  "ramp-response" 
imaging. 

In  1967  DeLorenzo  [Ref.  2]  presented  a  system  that  could  be  used  to  obtain  an  ap¬ 
proximate  measurement  of  the  impulse  response  of  scattering  bodies.  His  major  devel¬ 
opment  was  an  indoor  laboratory  scattering  range  which  facilitated  the  time  domain 
analysis.  He  showed  the  agreement  between  his  calculations  and  the  measurement  of  the 
impulse  response  of  a  rectangular  plate. 

Also  in  1967,  .Murchison  and  Falk  presented  a  hybrid  system  for  measurement  of 
impulse  responses  [Ref.  3).  This  measurement  system  was  separated  into  two  frequency 
ranges.  For  frequencies  below  4  Ghz,  they  used  an  impulse  generator  while  for  the  fre¬ 
quencies  above  4  Ghz,  they  used  a  swept  frequency  CW  source. 

In  1971  Baum  [Ref.  4]  demonstrated  theoretically  that  characteristic  resonances  can 
be  obtained  for  scattering  objects.  His  technique  was  called  the  Singularity  Expansion 
Method  (SEM).  It  was  postulated  that  the  resonances  could  be  obtained  from  the 
transient  scattering  measurement  for  a  complex  object. 

In  1973  Tesche  (Ref.  5]  applied  the  SEM  to  to  numerically  evaluate  complex  reso¬ 
nances  of  simple  metallic  scattering  objects. 
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In  I9ie  Blaricum  and  Mittra  [Ref,  6]  proposed  a  technique  for  extracting  the  poles 
and  residues  of  a  system  directly  from  its  transient  response.  Their  technique  provided 
a  number  of  numerical  advantages. 

In  1980  the  original  Transient  Electromagnetic  Scattering  Laboratory  (TESL)  be¬ 
came  operational  at  the  Naval  Postgraduate  School  (NPS).  Hammond  [Ref.  7]  devel¬ 
oped  this  laboratory  for  inverse  scattering  research  based  upon  synthesized  ramp 
responses.  This  laboratory  utilized  an  outside  ground-plane  range.  This  range  was 
limited  to  symmetric  scattering  objects  which  could  be  “mirror  imaged". 

In  1983  the  free-field  TESL  began  operation  at  NPS  [Ref  8].  This  laboratory  is  lo¬ 
cated  in  Spanagel  535,  The  TESL  facilitated  research  efforts  in  high  quality  transient 
scattering  measurements  (Ref  9]  used  for  radar  target  identification  and  characterization 
using  complex  pole  natural  resonances  [Ref.  10). 

Prior  to  the  current  undertaking,  McDaniel  [Ref  1!)  showed  good  agreement  be¬ 
tween  the  computation  and  the  measurement  of  the  impulse  responses  for  a  sphere  and 
a  thin  wire.  These  measurements  and  signal  processing  were  performed  using  the 
Tektronix  hardware  set  up,  as  originally  employed  in  the  1980  TESL.  In  1987  a  rev¬ 
olutionary  new  Hewlett  Packard  digital  processing  oscilloscope  became  available.  In 
addition,  computational  hardware  had  evolved  dramatically  since  the  early  1980  s.  This 
thesis  describes  the  implementation  of  this  new  technology  to  enhance  the  performance 
of  transient  scattering  and  signal  processing  at  the  NPS  TESL. 

B.  OVERVIEW  OF  THIS  EFFORT 

The  main  objective  of  this  thesis  was  to  incorporate  new  measurement  and  compu¬ 
tational  hardware  and  to  design  and  implement  new  software  to  automate  high-quality- 
transient  scattering  measurements.  In  this  thesis,  much  work  has  been  done  in  the  de¬ 
velopment  of  the  transient  scattering  work  station  which  has  the  ability  to  measure  and 
synthesize  the  scattering  impulse  response  of  scale  model  targets.  The  major  hardware 
includes  a  microcomputer,  a  digital  programable  oscilloscope  (DPO)  and  a  shielded 
anechoic  chamber.  Automatic  transient  scattering  measurements  have  been  imple¬ 
mented  by  utilizing  an  lEEE-bus  controller. 

With  the  new  HP54120T  DPO,  the  impulse  responses  of  scale  models  are  synthe¬ 
sized  at  a  high  signal-to-noise  ratio  (SNR).  The  acquisition  and  signal  processing  algo¬ 
rithms  run  much  faster  on  the  IBM  PC-AT  vis-a-vis  the  Tektronix  interpreted  BASIC 
computer  previously  used.  The  programs  now  use  QUICK  BASIC  and  FORTRAN 
compilers.  Validations  of  software  were  made  by  way  of  comparing  experimental 
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measurements  with  computations  for  simple  canonical  targets.  Various  scale  models 
were  measured  and  documented  in  detail.  Furthermore,  a  noise  calculation  was  provided 
for  further  research. 

This  thesis  is  broken  into  six  chapters.  Chapter  II  describes  system  facilities  and 
hardware.  This  includes  laboratory  description,  anechoic  chamber,  HP54120T  digital 
processing  oscilloscope  and  the  GURU  II  hardware  interface. 

Chapter  III  presents  the  theory  of  transient  scattering  measurements.  This  is  com¬ 
posed  of  the  free  field  range  system  representation,  mathematical  model  and  transient 
response  solution. 

Chapter  IV  deals  with  acquisition  algorithms  and  signal  processing.  This  chapter 
covers  operational  description,  GURU  II  software  interface  and  the  deconvolution  al¬ 
gorithm.  In  addition,  some  selected  outputs  are  illustrated. 

Chapter  V  emphasizes  the  experimental  results.  Two  types  of  targets  are  presented. 
One  is  a  canonical  target  and  the  other  is  a  scale  model  target.  The  noise  estimation  of 
the  system  is  also  discussed  in  this  chapter. 

Chapter  VI  summarizes  all  work  and  provides  some  future  considerations. 


II.  SYSTEM  FACILITIES  AND  HARDWARE 
A.  LABORATORY  DESCRIPTION 

Figure  1  on  page  5  shows  the  general  layout  of  the  scattering  range.  It  is  composed 
of 

•  Digital  Programmable  Oscilloscope  (DPO) 

•  Step  pulse  generator 

•  Low  loss  cables 

•  Broad  band  solid  state  power  amplifier 

•  Broad  band  transmitting  and  receiving  antennas 

•  Anechoic  chamber 

•  Personal  computer  (PC). 

Note  that  the  DPO  includes  the  sampling  and  pulse  generator  unit.  The  DPO  is 
used  as  both  the  master  pulse  source  and  as  the  receiver.  The  antenna  connection  cables 
were  obtained  from  the  DPO  accessory  kit  for  use  in  low  loss  network  measurements. 
They  are  flexible  for  hard-to-access  areas.  The  characteristics  are  shown  in  Table  21  and 
Table  22  on  page  78  in  Appendix  C. 

Transmitting  and  receiving  antennas  are  of  the  double-ridged  horn  type.  Their  fre¬ 
quency  range  is  1-12.4  GHz,  with  less  than  ldB  of  ripple  in  the  boresight  gain  over  this 
multi-octive  band.  The  Avantek  solid  state  power  amplifier  has  an  advertised  3dB 
bandwidth  of  2  to  6  GHz  but,  by  our  observations,  has  a  usable  bandwidth  of  1  to  7 
GHz.  A  15V  supply  is  needed  for  the  amplifier. 

The  anechoic  chamber  serves  two  purposes.  First,  it  is  used  to  isolate  the  target 
signal  from  the  electromagnetic  noise  and  interference  outside  of  the  chamber.  Secondly, 
the  absorbing  material  on  the  walls,  floor  and  ceiling  of  the  chamber  are  designed  to 
provide  minimal  reflections,  thus  approximating  unbounded  space.  The  computer, 
which  serves  both  acquisition  and  signal  processing  rolls,  is  an  1BM-AT  with  640  K 
RAM,  a  20  megabyte  hard  disk,  and  a  GPIB-PC2A  interface  card.  The  GPIB-PC2A 
was  added  to  the  PC  in  order  to  have  the  ability  to  talk  or  listen  to  the  HP54120T  DPO. 
The  PC  was  set  to  run  at  9  MHz.  It  has  two  floppy  disk  drives,  and  two  output  ports. 
One  is  a  parallel  port  for  the  printer  and  the  other  is  a  serial  port  for  the  plotter.  The 
color  graphic  board  is  NEC  GB1.  The  display  is  NEC  EGA  color  monitor. 
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B.  ANECHOIC  CHAMBER 

The  anechoic  chamber  is  located  in  room  535  of  Spanagel  Hall.  This  chamber  in¬ 
corporates  a  metallically  shielded  anechoic  enclosure  to  minimize  outside  EM  interfer¬ 
ence.  Chamber  dimensions  are  20  ft.  long  with  a  10  ft.  width  and  10  ft.  height.  Because 
targets  are  suspended  within  the  chamber,  there  is  no  requirement  that  they  have  mirror 
symmetry,  as  is  the  case  for  ground  plane  based  ranges.  In  the  measurements  that  were 
taken  for  this  thesis,  the  target  is  mounted  on  a  thin  styrofoam  column  at  a  distance  of 
8.5  feet  from  the  transmitting  and  receiving  antennas.  Both  receiving  and  transmitting 
antennas  are  placed  close  together  to  closely  approximate  monostatic  measurements. 
The  chamber  interior  is  covered  by  a  special  carbon  impregnated  foam  absorbing  mate¬ 
rial.  The  material  is  cut  into  longitudinal  wedges  along  the  side-walls,  floor  and  ceiling 
to  act  as  conduits  of  EM  energy  towards  the  back  wall.  The  back  wall  is  covered  by  18 
inch  long  pyramids  to  act  as  slow  absorbing  transitions  for  incident  fields.  Further  in¬ 
formation  about  this  free-field  TESL  can  be  found  from  Refs.  9,  11  and  12 

C.  HP54120T  DIGITAL  PROCESSING  OSCILLOSCOPE 

The  Hewlett  Packard  54120T  is  a  fully  programmable  20  GHz,  four  channel, 
digitizing  oscilloscope  with  a  nine  inch  color  display.  It  is  capable  of  automated  meas¬ 
urements,  digital  storage,  and  TDR  measurements.  This  DPO  has  only  been  in  pro¬ 
duction  since  mid- 1987  and  represents  the  state-of-the-art  in  sampling  technology. 
Previous  to  this  thesis  effort,  a  6  year  old  12  GHz  DPO  from  Tektronix  was  used  to 
perform  transient  scattering  measurements  in  the  TESL. 

The  HP-DPO  is  comprised  of  two  major  components.  These  are  the  HP54120A 
Digitizing  Oscilloscope  Mainframe  and  the  HP  541 21 A  Four  Channel  Test  Set.  It  also 
contains  extensive  self-tests  to  ensure  proper  functioning.  Both  HP  54120A  and  HP 
54121A  were  calibrated  together  as  a  system  called  the  HP  54120T.  The  DPO  can  be 
remotely  programmed  by  an  HP-IB  (IEEE-488)  controller  interface.  The  HP-IB  com¬ 
plies  with  the  IEEE  488.2  standard.  The  essential  characteristics  of  the  HP  541 20T  are 
summarized  in  Table  1  and  Table  2  on  page  7. 

It  should  be  noted  that  the  HP  541 21 A  4  Channel  Test  Set  is  sensitive  to 
electrostatic  discharge  (ESD).  Precautions  should  be  taken  to  keep  ESD  from  damaging 
the  equipment.  Details  can  be  found  in  Ref.  13. 

Waveforms  that  are  stored  in  the  DPO  memory  can  be  sent  to  either  an  HP  Thinkjet 
printer  or  an  HP  Plotter,  both  of  which  are  accessed  via  the  HP-IB  bus.  There  are  two 
ways  to  operate  the  instrument.  One  is  manual  operation  via  the  front  panel  and  the 
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Table  1.  HP54120T  CHANNEL  CHARACTERISTICS 


Bandwidth 

20/12.4  GHz 

Number  of  channels 

4 

Full  scale 

8  division 

Maximum  scale 

80  mV/division 

Minimum  scale 

1  mV/division 

Programable  dc  offset 

±  500  mV 

Input  maximum 

±  2  V  dc  4-  ac  peak 

Norminal  impedance 

50  ohm 

Connectors 

3.5  mm 

Percent  reflection 

less  than  5%  for  30  ps  risetime 

Table  2.  HP54120T  TDR  AND  TIMEBASE  CHARACTERISTICS 


Risetime 

less  than  45  psec  at  12.4  GHz  bandwidth  in  average 
mode 

Timebase  full-scale 

10  division 

Minimum  scale  factor 

10  psec/div 

Maximum  scale  factor 

Isec/div 

Time  interval  accuracy 

less  than  10  psec  +  0.1  %  of  reading 

other  is  to  operate  by  sending  HP-IB  bus  commands  via  the  computer.  A  major  efTort 
in  this  thesis  was  developing  programs  to  provide  easy  operator  control  and  signal  ac¬ 
quisition  of  the  DPO  using  the  computer.  These  algorithms  are  explained  in  detail  in 
chapter  IV. 

Ensemble  waveform  averaging  has  been  implemented  in  acquiring  scattering  signa¬ 
tures.  The  averaging  reduces  noise  by  l/JW ,  where  N  is  the  number  of  averages.  The 
minimum  discernible  signal  level  in  the  digital  processing  oscilloscope  system  is  approx¬ 
imately  35nV.  Typical  noise  amplitudes  for  the  DPO  sampling  unit  are  shown  in  Table 
3  on  page  8. 

Data  display  resolution  is  1024  points  horizontally  x  256  points  vertically  (for  dis¬ 
play  times  of  at  least  200  ps/div).  The  display  can  be  selected  in  persistence  or  average 


Table  3.  TYPICAL  NOISE _ _ 

20  GHz  bandwidth,  Avg  =  1 _  Noise  1.2  mV  (RMS) 

20  GHz  bandwidth,  Avg  =  256 _ Noise  (RMS) 

12.4  GHz  bandwidth,  Avg  =  1  Noise  500/*F( RMS) 

12.4  GHz  bandwidth,  Avg  =  256 _ Noise  35m F (RMS) 

12.4  GHz  bandwidth,  persistence  Noise  400mF  (RMS) 

mode  and  the  number  of  averages  can  be  specified  as  a  power  of  2,  up  to  2048.  The  data 
output  transfer  rate  and  the  data  record  lengths  are  shown  in  the  following  table. 

Table  4.  DATA  OUTPUT  TRANSFER  RATE  AND  DATA  RECORD  LENGTHS 


10  psec/div  <  time  div  <  20  psec/div 

100,  or  400  points  record 

20  psec'div  <,  time  div  <  50  psec/div 

100,  400  or  800  points  record 

50  psec  div  <.  time  div  <  200  psec/div 

100,  500  or  1000  points'record 

200  psec/div  <,  time/div  <  1  sec,  div 

128,  256,  500,  512  or  1024  points/record 

Data  output  transfer  rate 

115  kbytes/s 

For  the  time  domain  reflectometer  (TDR)  system,  the  risetime  is  an  adjustable  value 
based  upon  the  timebase  setting.  The  minimum  is  10  ps  or  .08x  time/div,  whichever  is 
greater.  The  maximum  is  5  x  time/div.  Usually  channel  1  is  used  for  reflection  meas¬ 
urements,  however,  in  this  thesis,  it  was  used  as  the  source  generator. 

The  digitizer  uses  a  12-bit  successive  approximation  A/D  converter.  The  full-scale 
range  of  the  A  D  is  640  mV.  The  least  significant  bit  (LSB)  of  the  A/D  converter  equals 
250  txV.  This  gives  one  part  in  2560,  or  more  than  11  bits  of  resolution.  Averaging  can 
stretch  out  the  resolution  to  32  nV.  This  improves  the  resolution  to  around  14  bits. 

The  signal  is  sampled  and  digitized  at  a  rate  determined  by  the  trigger  rate,  repe¬ 
tition  rate,  time  base  range,  display  mode,  and  number  of  channels  turned  on.  Further 
information  concerning  the  characteristics  of  the  HP  54120T  can  be  found  from  Ref. 
14. 

D.  GURU  II  HARDWARE 

The  Textronix  supplied  General-purpose  interface  bus  User's  Resource  Utility 
(GURU  II)  is  a  combination  hardware  and  software  package  that  allows  the  use  of  any 


of  the  IBM  personal  computers  (PC),  and  some  PC  compatibles,  to  control  GPIB 
programable  instruments.  In  this  section,  the  focus  is  on  the  hardware.  The  software 
is  explored  in  the  next  chapter. 

The  hardware  consists  of  the  PC2A  GPIB  interface  card  that  plugs  into  the  PC  main 
board,  and  the  GPIB  cable.  The  GPIB  PC2A  interface  card  is  an  8-bit  half-size  board 
which  performs  as  a  talker,  listener  and  controller.  It  also  can  be  implemented  with  the 
full  range  of  talker,  listener,  serial  poll,  service  request,  and  remote  programming  func¬ 
tions  for  as  many  as  14  devices  on  the  bus.  It  can  perform  as  a  complete  controller. 
Using  the  PC2A,  data  transfer  between  the  system  memory  and  the  GPIB  are  performed 
at  rates  of  more  than  300  kilobytes  per  second.  The  PC2A  is  compatible  with  all  revision 
levels  of  the  IEEE-488  standard  including  the  HP-IB. 

The  GPIB-PC2A  interface  card  consists  of  these  major  sections: 

•  address  decoding 

•  buffering  and  data  routing 

•  interrupt  arbitration 

•  direct  memory  access  (DMA)  arbitation. 

•  configuration  switches  and  jumpers 

•  GPIB-adapter-TLC  (talker,  listener/controller) 

•  time  of  day  clock  with  battery  back-up. 

The  address  decoding  monitors  the  address  lines  to  recognize  when  the  GPIB-PC 
I/O  address  is  present  on  the  computer  I/O  channel  and  enables  read  and  write  access 
to  the  GPIB  adapter  register.  The  buffering  and  data  routing  handles  data  transfer  be¬ 
tween  the  IBM  PC-AT  I/O  channel  and  the  GPIB  adapter  through  a  bidirectional 
internal  data  bus. 

Interrupt  arbitation  recognizes  when  interrupts  have  been  enabled  or  disabled  and 
passes  or  inhibits  them  accordingly.  DMA  arbitation  recognizes  when  DMA  operations 
are  enabled  or  disabled,  and  when  the  last  transfer  has  taken  place.  It  also  routes  the 
DMA  request  and  acknowledges  signals  to  the  selected  DMA  channel.  The  DMA  is 
used  to  transfer  data  directly  between  a  peripheral  and  memory  without  CPU. 

The  IEEE-bus  controller  forms  an  interface  between  the  computer  and  the 
IEEE-488  bus.  Conceptually,  the  IEEE-488  bus  behaves  in  a  way  very  similar  to  the 
system  bus.  Hewlett-Packard  originally  devised  it  (Ref.  15]  and  intended  it  to  link  to¬ 
gether  programmable  instruments  in  a  laboratory  or  industrial  environment.  By  con¬ 
trolling  test  equipment  and  measuring  devices  from  the  IEEE-488  bus,  implementing  an 


automatic  testing  station  is  possible.  A  system  under  test  is  connected  to  the  test 
equipment  and  measuring  devices.  The  computer  configures  all  equipment  via  the 
IEEE-488  bus  and  then  reads  the  test  results  from  the  same  bus.  The  GP1B  has  the 
following  limitations: 

•  Half  of  the  devices  attached  to  the  GPIB  must  have  power  on. 

•  The  total  length  of  the  bus  cable  must  not  exceed  20  Meters. 

•  No  more  than  15  devices  may  be  attached  to  one  bus. 

More  detailed  information  can  be  found  in  Ref.  16. 


III.  THEORY  OF  TRANSIENT  SCATTERING  MEASUREMENTS 

A.  FREE-FIELD  RANGE  SYSTEM  REPRESENTATION 

The  physical  model  of  the  system  free  field  scattering  range  was  first  presented  by 
Morgan  in  Ref.  12.  The  system  block  diagram  is  showm  in  Figure  2  on  page  12.  Im¬ 
portant  components  and  parameters  of  the  system  are  summarized  in  the  frequency  do¬ 
main  representation  in  the  following  table. 


Table  5.  FREQUENCY  DOMAIN  FUNCTIONS 


Pulse  generator,  driving  the  transmitting  antenna 

Xif) 

Transmitting  antenna  response 

HrV) 

Receiving  antenna  response 

HJJ) 

The  multiple  scattering  between  target  and  chamber  wall 

E3SB9ZH 

Direct  couple  of  receiving  and  transmitting  antenna 

HdV) 

The  multiple  scattering  between  calibration  sphere  and 
absorber 

him 

The  multiple  scattering  between  target  and  absorber 

HIM 

The  calibration  sphere  transfer  function 

mv) 

The  target  transfer  function 

HIV) 

Thermal  noise 

m 

The  transient  pulse  source  is  the  step  generator  in  the  TDR  system.  This  generator 
not  only  supplies  the  low-level  transmitted  pulse,  before  amplification,  but  also  acts  to 
coherently  trigger  the  sampling  circuits  of  the  DPO  receiver.  It  thus  is  akin  to  the 
“master  oscillator"  in  a  coherent  pulsed  CW  radar  system.  The  risetime  of  the  step 
generator  is  approximately  25  psec.  This  step  pulse  is  first  amplified  and  filtered  by  a 
broad-band  GaAs  FET  amplifier. 

The  amplified  pulse,  which  no  longer  resembles  a  step  waveform,  is  then  fed  through 
the  transmitting  antenna  which  radiates  into  the  anechoic  chamber.  The  radiated  field 
emanating  from  the  antenna  is  approximately  an  angular  weighted  spherical  wave.  At 
the  target  position  the  incident  field  is  a  quasi-plane  wave.  Because  of  the  close  prox¬ 
imity  of  the  transmitting  and  receiving  antennas,  some  strong  coupling  of  energy  is  in¬ 
duced  in  the  receiving  antenna. 
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This  undesirable  signal  can  be  seen  clearly  if  the  time  delay  setting  on  the  DP O  is  ad¬ 
justed  to  begin  many  nanoseconds  before  the  target  response  time-window.  Therefore, 
the  time  setting  delay  must  be  checked  when  the  target  is  present. 

The  initial  calibration  of  the  beginning  of  the  target  time-window  was  performed  by 
using  a  rectangular  copper  plate,  1  ft  wide  and  1.5  ft  high,  as  the  test  scatterer.  The  back 
scattering  was  observed  with  the  copper  plate  turned  broadside  to  the  centroid  of  the 
horn  positions  and  placed  where  the  target  would  be  mounted.  The  copper  plate  scat¬ 
tering  signature  was  so  large  that  it  could  be  clearly  discerned.  The  results  of  this 
measurement  indicated  a  DPO  time  delay  setting  of  38  nsec  was  needed. 

Multiple  scattering  between  the  target  and  the  chamber  absorber  is  represented  by 
the  two  way  interaction  arrow'  connecting  Hdf)  to  Hs[f ).  The  total  noise  of  the  system 
A ’(/),  is  due  to  both  thermal  emission  from  the  chamber  and  receiver  noise  in  the  sam¬ 
pling  front  end  of  the  DPO.  These  sources  are  combined  together  and  denoted  as  A'„(/) 
for  the  n-th  measurement  that  is  made.  To  understand  the  theory  of  transient  scattering 
measurements  adequately,  the  mathematical  model  will  be  presented  in  the  next  section. 

B.  MATHEMATICAL  MODEL 

The  five  key  parameters  in  Table  6  on  page  14  play  a  major  role  in  both  the  physical 
meaning  and  mathematical  sense  of  the  acquisition  and  signal  processing  algorithms. 
These  five  parameters  can  be  related  to  the  system  representation,  described  in  the  pre¬ 
vious  chapter,  by  the  following  three  equations: 

= im  n*(/)  m  [ihV) + HcV) + HsV) + /&(/)] + m  pi) 

v2(f)  -  im  im  xn)  + //c(/)j  +  a2oo  (3.2) 

y3(f>  -  im  im  m  [im + "d/) + "Jw + "Id/)] + m-  (3.3) 

These  three  equations  represent  the  three  respective  measurements  that  are  neces¬ 
sary  for  processing  the  scattering  signature  of  a  given  target.  As  shown  below,  these 
measurements  are  for  a  calibration  sphere,  whose  measured  scattering  can  be  compared 
to  an  accurate  calculation;  the  "background"  which  denotes  no  target  present,  and  the 
actual  target  scattering  measurement. 

Equation  3.1  states  that  the  response  to  the  amplified  pulse  generator,  X(f).  due  to 
the  calibration  sphere  and  the  chamber  is  equal  to  the  product  of  the  pulse  excitation, 
transmitting  antenna  transfer  function  //,(/),  receiving  antenna  transfer  function  //*(/) 
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Table  6.  MATHEMATICAL  NOTATION 


The  calibration  measurement  waveform 

»',(/) 

The  background  measurement  waveform 

YJJ) 

The  target  measurement  waveform 

YJJ) 

The  subtracted  calibration  waveform 

YJJ) 

The  subtracted  target  waveform 

YJJ) 

times  the  sum  of  the  antenna  coupling  HA(J)  plus  the  calibration  sphere  7/f(/)  and  the 
multipath  multiple  scattering  Added  to  this  product  is  the  thermal  noise  A„(/). 

Equation  3.2  represents  the  background  measurement,  with  no  target  present.  The 
system  response  is  similar  to  that  previously  considered,  but  with  the  absence  of  the 
transfer  function  for  the  calibration  sphere.  Equation  3.3  expresses  the  target  measure¬ 
ment  in  the  frequency  domain.  The  intepretation  follows  the  same  pattern  as  equation 
3. 1  except  the  target  transfer  function,  ///(/),  appears  instead  of  that  for  the  calibration 
sphere.  The  multipath  multiple  scattering  in  equation  3.3  also  changes  and  is  denoted 

by  irsc(f). 

The  three  equations  as  stated  above  are  three  crucial  measurements  in  the  transient 
scattering  problem.  I  he  backscatter  measurement  with  the  canonical  (calibration)  ob¬ 
ject  }',(/),  measurements  without  target  (background)  Y,{J)  and  measurements  with  the 
desired  target  1  }{f)  are  obtained  directly  using  the  time  domain  acquisition  program  in 
Appendix  B.  This  program  was  written  in  BASIC  as  part  of  this  thesis  effort.  To  speed 
up  the  execution  of  the  program  by  a  factor  of  about  fifty  times  vis-a-vis  direct  inter¬ 
pretation,  it  has  been  compiled  using  the  Microsoft  Quick  Basic  Compiler. 

After  acquiring  accurate  data,  the  next  step  is  to  process  the  data  and  extract  the 
impulse  response.  This  is  the  topic  of  the  next  section. 

C.  TRANSIENT  RESPONSE  SOLUTION 

The  three  measurements  that  were  just  considered  form  the  basis  for  estimating  the 
“smoothed"  impulse  response  of  the  target.  The  band-limited  nature  of  the  measure¬ 
ments  does  not  allow  synthesis  of  the  true  infinite  bandwidth  impulse  response  which  is 
due  to  an  incident  delta  function  plane  wave.  Instead,  the  "smoothed"  term  means  that 
we  can  estimate  the  transient  target  response  due  to  some  predefined  incident  signal  (a 
double-Gaussian  shaped  waveform  is  used  here)  whose  significant  spectral  content  docs 
not  exceed  the  bandwidth  of  the  measurements.  Aside  from  the  currently  used  solid- 
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state  power  amplifier,  the  bandpass  range  of  the  DPO-antenna  system  is  about  1  GMz 
to  12  GHz.  The  spectral  bottleneck  is  due  to  the  GaAS  FET  amplifier,  whose  passband 
is  about  1  GHz  to  7  GHz. 

In  order  to  synthesize  the  smoothed  transfer  function  of  the  target  ///(/),  the  first 
step  is  to  subtract  the  "clutter"  portion  of  the  responses  in  equations  3.1  and  3.3  which 
is  embodied  in  equation  3.2.  This  subtraction  results  in  equations  3.4  and  3.5  below: 


s 

R 

S 

1 

s 

(3.4  a) 

Y,{f)  -  Hjifi  Hr{J)  X{j)  +  /&(/)]  +  A'4(/) 

(3.4  b) 

> 

*» 

s 

II 

s 

1 

> 

s 

(3.4c) 

i's(/)-*3(/;->j(/) 

(3.5a) 

y5{/) = im  //*(/)  x(j)  [nim + y/jc(/)] + Aj(/> 

(3.56) 

(3.5c) 

Equations  3.4  and  3.5  are  the  very  critical  steps  in  obtaining  the  smoothed  transient 
response  of  the  target.  The  critical  nature  of  the  subtractions  has  to  do  with  coherency 
in  time  of  the  subtracted  clutter.  The  calibration  sphere  measurement  }',(/)  had  the 
background  measurement  1  \[J)  subtracted  in  order  to  eliminate  the  clutter.  The  clutter 
signal  is  that  which  is  present  aside  from  the  target  and  is  composed  of  two  main  con¬ 
stituents.  One  is  the  direct  coupling  of  the  receiving  and  transmitting  antennas,  //*(/), 
and  the  other  is  direct  scattering  (as  opposed  to  multipath)  from  the  absorber  in  the 
chamber,  lic (/).  Since  the  measurements  are  made  directly  in  the  time  domain,  the  sub¬ 
tractions  are  done  in  time  domain.  In  a  similar  fashion,  equation  3.5  can  be  obtained. 

The  results,  )'„(/)  and  15(f),  are  then  transformed  into  frequency  domain  by  discrete 
fast  Fourier  transform  (  FFT  ).  This  leads  to  the  next  step,  which  is  to  form  an  optimal 
deconvolution  estimator  for  target  response  given  by 


A'oOO  //[(/)  - 


m  Km  ly  e/) 

l>'«(/)l2  +  C 


(3.6) 


where  >  *(/)  is  the  complex  conjugate  of  }'„(/),  YSC(J)  is  the  computed  calibration  sphere 
and  C  is  the  smoothing  parameter. 
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YscU)  =  Xdj)HCs(f) 

Xo(0  *-  FFT  -*  XcfJ) 


(3.7) 


(3.8) 

j t(k)  -  c,  <r 3,1  (/_  °2  +  c2 '°)J,  (3.9) 

where,  for  a  specified  a,  and  otj  in  this  "Double-Gaussian"  (  DG  )  pulse,  the  C,  and  C, 
are  chosen  such  that 

x0(t0)  =  1  (3.10) 

and 


/•  +00 

j_^Jr0(/y/  =  °.  (3.11) 

Each  of  the  two  Gaussian  waveform  in  equation  3.9  will  be  specified  here  by  a  "10% 
pulse  width",  defined  by  A i,  such  that 

e-*>—r~  =  0.1,  (3.12a) 


which  gives 


ak 


4 

(A  fir)2 


ln(10) 


9.21 

(A/*r  ' 


(3.126) 


The  estimator  in  equation  3.6  is  due  to  Riad  [Ref.  17]  and  results  from  an  optimal 
least-square  fidelity  criterion  for  deconvolved  target  response.  The  result  of  equation 
3.6,  Xdf)  Hi ,  represents  the  frequency  domain  scattering  response  of  the  measured  tar¬ 
get  to  an  incident  plane  wave  whose  Fourier  transform  is  X0{f).  This  incident  waveform 
is  specified  by  numerically  computing  the  scattering  response  of  the  calibration  sphere 
to  it.  This  numerically  generated  response  of  the  calibration  sphere  is  embodied  in 
Yscif)  >  as  shown  in  equation  3.7.  It  is  obtained  from  the  SPRSCT  program  which  uses 
the  magnitude  and  phase  data  from  the  MIE  series  program.  Both  of  these  programs 
were  written  by  Morgan. 

Both  subtraction  and  optimal  deconvolution  algorithms  are  implemented  in  the 
Deconvolution  Program  in  Appendix  A.  The  procedure  to  operate  this  program  will  be 
explained  in  chapter  IV. 
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IV.  ACQUISITION  AND  SIGNAL  PROCESSING  ALGORITHMS 


A.  GURU  II  INTERFACE 
1.  GPIB  Principal 

This  section  emphasizes  the  GURU  II  software  interface.  In  order  to  under¬ 
stand  the  Acquisition  program  operation,  the  GPIB  principle  has  to  be  understood.  The 
General  Purpose  Interface  Bus  uses  a  twenty-four  conductor  cable.  Sixteen  of  these  are 
used  to  send  the  signal.  Data  is  sent  back  and  fourth  between  the  listener  and  the  talker 
in  a  bi-directional  fashion.  This  means  that  both  transmitted  and  received  data  travel 
on  the  same  line.  Eight  lines  are  used  for  8  bits  of  parallel  data  transfer.  A  byte  is  sent 
by  serial  mode  but  each  bit  is  sent  in  parallel  mode.  Three  lines  act  as  the  coodinator 
between  the  talker  and  the  listener.  These  lines  are  assigned  their  duty  as  follows: 

•  no  data  available  line 

•  data  not  ready  line 

•  data  not  accepted  line. 

The  talker  on  the  bus  controls  the  data  available  on  the  line  and  waits  for  the 
listener  s  response  on  the  "not  ready  for  data"  line  or  "data  not  accepted"  line.  The  re¬ 
maining  five  lines  are  the  attention  line,  the  interface  clear  line,  the  remote  enable  line, 
the  service  request  line  and  the  end  or  identify  line.  When  the  attention  line  is  active, 
the  talker  and  the  listener  will  be  specified.  The  controller  has  full  access  to  five  of  these 
lines  and  has  the  capability  to  directly  address  other  devices  to  command  them  to  talk 
or  listen.  The  interface  clear  line  is  used  to  set  the  turn-on  condition.  When  the  remote 
enable  line  is  asserted,  all  devices  on  the  bus  are  under  computer  control  and  nothing 
external  can  interfere.  The  service  request  line  is  used  by  a  device  to  indicate  a  require¬ 
ment  for  attention  and  to  request  an  interruption  of  the  current  sequence  of  events. 
There  are  two  methods  for  the  controller  to  determine  the  device  and  the  message.  One 
is  the  serial  poll  while  the  other  is  the  parallel  poll.  The  polls  obtain  numerical  codes 
which  show  a  specific  event.  The  serial  poll  is  a  mechanism  that  the  controller  uses  to 
obtain  status  from  individual  devices.  As  each  device  is  polled,  it  returns  a  status  byte 
which  controller  decodes  to  determine  the  device  requesting  service.  The  parallel  poll  is 
a  mechanism  for  accepting  and  decoding  the  general  status  of  as  many  as  eight  devices 


simultaneously.  The  end  or  identify  line  is  used  by  the  talker  to  indicate  the  end  of  the 
message. 

2.  Basic  Software  Requirements 

There  are  3  files  which  the  GURU  requires  in  order  to  run  the  application  pro¬ 
gram.  These  files  are 

•  GPIB.COM 

•  IBCONF.EXE 

•  CBIB.OBJ. 

The  primary  GPIB  device  drivers  are  read  in  via  the  CONFIG.SYS  file  at  boot-up  time. 
This  file  is  in  the  root  director}'  and  must  include  the  following  DOS  statement: 

•  DEVICE  =  GPIB.COM. 

This  causes  the  GPIB  driver  in  GPIB.COM  to  be  loaded  when  the  system  is  booted. 

The  IBCONF.EXE  is  the  executable  program  to  configure  the  GPIB.COM  file 
containing  the  device  drivers  for  each  instrument.  HPDEV  is  the  name  of  the  device 
driver  for  all  Hewlet-Packard  instruments,  including  the  HP  54120T,  whose  decimal  pri¬ 
mary  address  is  17.  The  IBCONF.EXE  (Interface  Bus  Configure)  program  allows  a 
match  between  the  softwere  and  the  hardware.  The  program  prompts  with  menus  of 
details  about  the  hardware  parameters  associated  with  the  interface  card.  This  device 
for  Hewlett-Packard  is  already  configured  in  the  system.  When  the  IBCONF.EXE  is 
run,  the  device  drivers  are  mapped  for  each  instrument  in  the  display.  By  selecting  the 
HPDEV,  and  editing  it,  the  characteristics  of  the  HPDEV  show  on  the  screen.  This  al¬ 
lows  the  operator  to  change  the  characteristics  to  be  compatible  with  the  hardware. 
After  reconfiguration  of  the  GPIB.COM  file,  the  operator  must  return  to  the  operating 
system  and  then  re-boot  the  PC  for  the  new  parameters  to  be  in  effect. 

CBIB.OBJ  is  the  object  file  which  is  linked  with  any  compiled  BASIC  program 
which  needs  to  interface  with  the  GPIB.COM  device  drivers.  In  our  case,  the  Acquisi¬ 
tion  Program,  written  in  GW-BASIC,  is  first  compiled  using  the  Microsoft  QUICK 
BASIC  compiler.  This  process  creates  an  object  (.OBJ)  file.  The  CBIB.OBJ  object  file 
is  then  linked  with  the  Acquisition  object  file.  The  final  result  is  an  executable  (.EXE) 
file. 

B.  ACQUISITION  ALGORITHM 

The  Acquisition  program  was  written  in  Microsoft  GW-BASIC.  The  hardware  is 
interfaced  with  this  program  by  way  of  the  CBIB  object  file.  This  file  comes  with  the 


18 


GURU  II  software.  After  compiling  the  source  code,  the  object  file  of  the  Acquisition 
has  to  be  linked  to  the  CBIB  object  file  in  order  to  complete  the  creation  of  the  Acqui¬ 
sition  executable  file.  To  run  the  program,  the  Brun3087  execution  file  is  needed  in  the 
same  DOS  Directory.  This  file  is  supplied  with  Microsoft  QUICK  BASIC.  The  Ac¬ 
quisition  program  uses  a  friendly  and  interactive  menu  format.  The  program  begins  with 
a  displayed  welcome  message,  including  the  time  and  date.  This  screen  also  provides  a 
brief  procedure  on  how  to  measure  the  transient  scattering  response.  The  screen  gives 
operator  the  choice  of  using  the  system  default  or  the  last  response.  The  system  default 
was  designed  for  first  time  operators.  The  last  response  was  designed  for  experienced 
operators.  A  main  menu  is  displayed  on  the  screen  after  the  operator  decides  whether 
to  use  the  system  default  or  the  last  response. 

The  main  menu  contains  many  importance  parameters.  A  first  item  is  the  number 
of  data  points  in  the  acquired  time  series.  This  number  has  to  be  a  power  of  2,  but  not 
more  than  1024.  This  limitation  could  be  modified  in  the  program  by  defining  a  larger 
array  size  up  to  2048  points.  The  second  item  is  the  number  of  sub-averages.  This  pa¬ 
rameter  determines  the  number  of  values  to  average  for  each  time  point  when  the 
HP54120T  is  in  the  average  mode.  This  parameter  has  to  be  an  integer  from  1  to  2048, 
also  in  a  power  of  2.  The  third  item  is  the  number  of  data  blocks.  This  determines  the 
number  of  sets  of  data  in  one  record.  The  fourth  and  the  fifth  items  are  a  user  supplier 
identification  for  the  target  and  the  date  of  measurement.  The  sixth  item  tells  the  oper¬ 
ator  about  the  typical  sampling  time  interval.  The  seventh  parameter  is  the  type  of  the 
waveform  being  acquired  on  the  measurement  being  made.  As  discussed  in  the  previous 
chapter,  there  are  three  scattering  waveforms  which  need  to  be  acquired  for  each  target 
being  considered:  calibration  (sphere),  background  and  actual  target.  The  eighth  item 
tells  the  operator  the  time  window  in  nanoseconds.  The  ninth,  tenth  and  eleventh  items 
are  the  maximum  vertical  setting,  the  data  file  name  output  and  the  auto  time  setting. 
At  the  bottom  of  the  menu,  the  operator  is  asked  if  he  wants  to  change  the  items. 

The  second  menu  appears  on  the  screen  when  the  operator  decides  to  change  some 
item  on  the  main  menu.  This  "overlay"  menu  contains  all  changeable  parameters  and 
allows  a  change  to  these  parameters  interactively.  When  the  changes  are  completed,  the 
operator  can  go  back  to  the  main  menu  by  typing  the  letter  "R".  This  brings  the  oper¬ 
ator  to  the  main  menu  and  displays  the  new  item  on  the  menu.  When  everything  is 
ready,  the  operator  selects  "no  change"  to  begin  to  automate  the  DPO.  The  program 
starts  by  initializing  the  DPO  and  setting  all  specified  parameters  as  directed  by  the  main 


menu.  The  program  stops  when  the  DPO  has  finished  the  setting  operation.  The  pro¬ 
gram  then  asks  for  the  bandwidth  of  the  measurement.  There  are  two  choices:  20  Ghz 
or  12.4  Ghz.  To  protect  from  errors  in  the  initial  setting,  the  program  stops  for  the  final 
check.  As  soon  as  the  operator  initiates  the  acqisition  cycle,  the  beginning  time  shows 
on  the  screen.  For  a  typical  acquisition  composed  of  2048  averages  and  1  data  block, 
it  takes  about  7  minutes  to  finish  the  operation.  In  the  meantime,  the  computer  waits 
for  the  DPO  to  provide  the  raw  data.  When  the  DPO  finishes  this  operation,  the  ending 
time  is  displayed  and  the  PC  opens  a  temporary  hard-disk  file  for  storing  the  raw  data. 
Before  the  operation  is  ended,  the  header  which  is  to  be  placed  on  the  data  file  is  dis¬ 
played  on  the  screen.  The  operator  can  then  observe  any  error  in  the  data  set,  should 
there  be  any. 

The  header  is  transferred  as  a  separate  file.  The  header  provides  important  charac¬ 
teristics  of  the  data,  such  as  format  type,  number  of  points,  number  of  averages,  x  in¬ 
crement,  x  origin,  x  reference,  y  increment,  y  origin,  y  reference  and  range  setting.  All 
of  these  parameters  allow  the  program  to  recover  meaningful  data  from  the  raw  data. 
The  remainder  of  the  program  constructs  the  standard  data  file  structure  and  the  hard 
copy  output. 

At  the  end  of  the  program,  there  is  a  pause  with  a  query  regarding  the  operator's 
intention  to  continue  with  an  additional  measurement.  To  complete  the  acquisition  cy¬ 
cle  for  one  target,  the  three  measurements  previously  discussed  have  to  be  made  in  order. 
These  three  files  will  then  be  supplied  to  the  deconvolution  program  in  the  next  proc¬ 
essing  step.  For  acquiring  more  than  three  measurements,  the  background  must  be 
measured  for  each  target  so  it  can  be  subtracted  from  the  respective  target  measurement. 
The  background  measurement  occurs  either  immediately  before  or  after  the  target 
measurement.  This  reduces  errors  induced  by  long  term  nonstationary  drifts  in  the  sys¬ 
tem.  The  need  for  such  temporal  coherence  of  the  target-background  measurements  has 
been  verified  experimentally  during  the  course  of  this  work.  Typical  outputs  are  as 
shown  in  Figure  3  on  page  21. 

C.  DECONVOLUTION  ALGORITHM 

The  deconvolution  program  was  written  in  Ryan-McFarland  (RM)  FORTRAN  in 
order  to  improve  the  run  time.  The  source  code  is  standard  ANSI  77  FORTRAN  so  it 
may  be  recompiled  under  virtually  any  PC  FORTRAN  compiler  such  as  Microsoft  or 
Lahey.  This  program  was  originally  written  in  Tektronix  BASIC  by  Morgan  [Ref  11] 
and  was  restricted  to  use  on  a  Tektronix  4052  microcomputer. 
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Figure  3.  Deconvolution  input  waveform 


The  use  of  a  standard  FORTRAN-77  compiler  generalizes  the  program  so  it  can  be 
run  on  any  IBM  PC  compatible.  This  algorithm  follows  the  mathematic  deconvolution 
as  described  in  chapter  III.  The  program  has  graphic  capabilities  which  are  supplied  by 
GRAFMATIC  software.  The  procedure  to  compile  the  program  is  described  in  the  RM 
FORTRAN  manual.  This  program  must  be  linked  to  the  GRAFMATIC  library  in  or¬ 
der  to  run  properly. 

The  program  begins  by  displaying  the  standard  input  data  file.  This  data  file  format 
has  been  standardized  for  use  with  all  of  the  software  support  systems.  The  header  of 
each  data  file  contains  the  important  information  such  as  identification,  type  of  meas- 
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urement,  date  of  measurement,  time  of  experiment,  the  maximum  time  scale,  data  block 
and  number  of  points.  Types  of  measurement  are  as  explained  in  the  previous  section. 
The  first  data  file  needed  for  the  deconvolution  program  is  that  for  the  computed  scat¬ 
tering  by  the  calibration  sphere  target.  This  computation  data  is  obtained  from  the  MIC 
series  program  and  SPRCT  program,  written  by  Morgan,  as  described  in  the  previous 
chapter.  The  input  data  required  for  these  programs  are  summarized  in  Table  7  on  page 
23  . 

The  sphere  scattering  computation  is  transformed  into  the  frequency  domain  and 
displayed  on  the  screen  in  order  to  observe  the  bandwidth.  The  bandwidth  of  the  com¬ 
putation  has  to  be  approximately  the  same  as  the  estimated  bandwidth  of  the  measure¬ 
ment.  If  the  bandwidth  is  not  appropriate,  the  operator  should  to  terminate  the 
program  via  a  Ctl-C  and  run  the  MIC  series  followed  by  SPRCT  program  to  provide 
appropriate  data. 

In  the  next  step,  the  program  asks  for  the  calibration  and  background  file  names. 
Both  data  are  read  in  and  plotted  on  the  same  axis  in  order  to  observe  if  time  shifting 
has  occurred,  and  if  so,  to  be  able  to  estimate  how  much.  The  operator  is  then  allowed 
to  shift  one  waveform  cither  left  or  right  and  then  re-observe  the  overlayed  waveforms. 
Once  a  satisfactory  alignment  is  observed,  the  calibration  waveform  is  subtracted  from 
the  background  waveform  and  the  difference  waveform  is  transformed  into  the  frequency 
domain. 

The  program  then  stops  and  asks  for  the  target  file  name.  It  then  asks  if  the  oper¬ 
ator  desires  to  use  the  same  background  or  another,  depending  on  which  background 
was  closer  to  the  target  completion  acquiring  time.  In  a  similar  fashion,  the  target  data 
and  the  background  data  are  plotted  on  the  same  axis.  If  there  is  an  observed  drift,  the 
program  allows  the  operator  to  shift  the  target  waveform  to  the  left  or  to  the  right  until 
both  target  and  background  arc  lined  up.  Then  the  subtraction  can  be  performed.  The 
difference  waveform  is  then  transformed  into  the  frequency  domain.  Typical  scattering 
data  being  processed  by  deconvolution  program  is  shown  in  Figure  3  on  page  21.  The 
overlay  and  shift  is  shown  in  Figure  4  on  page  23  and  the  subtracted  waveform  is  shown 
in  Figure  5  on  page  24.  The  target  for  this  example  was  a  10  centimeter  long  "thick" 
copper  wire,  and  the  calibration  target  was  a  standard  0.1025  meter  radius  aluminum 
sphere.  After  the  subtraction,  the  difference  calibration  waveform  and  the  difference 
target  waveform  were  tranformed  into  the  frequency  domain. 
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Figure  4.  Overlay 


Table  7.  CALIBRATION  SPHERE  PARAMETERS 


Sphere  radius  in  meters 

0.1025 

Distance  from  sphere  in  meters 

2.5 

Bistatic  angle  in  degrees 

3.0 

Scattering  plane 

L  plane 

Time  window  in  nanoseconds 

20.0 

Nyquist  frequency  in  Ghz 

12.775 

No.  frequency  points 

257 

23 


SUBTRACTED  WAVEFORM 


Figure  5.  Subtraction  waveform 

The  three  frequency  domain  waveforms  were  formed  into  Riad's  optimal  deconvo¬ 
lution  estimator  given  in  equation  3.6.  The  smoothing  constant  "C"  is  selected  by  the 
operator.  This  constant  is  set  to  optimize  the  denominator  threshold  level  in  the  fre¬ 
quency  domain.  If  it  is  set  too  large,  oversmoothing  will  result  which  reduces  the  spec¬ 
tral  content  in  both  the  low  and  high  frequency  regimes.  If  the  threshold  is  set  to  zero, 
a  "naive"  deconvolution  results.  This  provides  vastly  amplified  noise  in  the  signal  and 
typically  looks  like  "garbage”.  A  typical  value  for  C  is  0.1.  Riad's  method  yields  an 
estimator  for  the  scattered  w'aveform  of  the  target  due  to  an  incident  plane  wave  signal 
having  signal  shape  equal  to  that  which  illuminates  the  computed  calibation  sphere  us¬ 
ing  program  SPRCT.  This  frequency  domain  estimator  is  transformed  to  the  time  do- 
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main  and  the  result  is  displayed.  At  the  end  of  the  program,  the  operator  may  try 
another  smoothing  parameter.  If  the  result  satisfies  the  operator,  it  can  be  stored  in  the 
file  whose  file  name  was  provided  by  the  operator.  The  program  also  allows  the  proc¬ 
essing  of  another  target  directly  without  rereading  the  computed  data  and  the  calibration 
data. 


D.  EXAMPLE  OUTPUT 

The  output  of  the  deconvolution  program  is  in  both  the  time  domain  and  frequency 
domains,  but  the  following  example  shows  only  the  time  domain  output.  This  output 
can  be  stored  on  a  hard  disk  or  a  floppy  disk. 


Figure  6.  Thick  wire  measurement 
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The  impulse  response  shown  in  Figure  6  belonged  to  the  signature  of  the  thick  copper 
wire  having  a  90  degree  incident  wave  and  a  90  degree  reflected  wave.  The  number  of 
averages  used  in  the  measurement  was  2048.  The  smoothing  was  C  =  0.1.  There  was 
no  time  drift  observed.  One  data  block  was  used.  Other  parameters  can  be  found  in  the 
following  table. 


Table  8.  90  DEGREE  THICK  WIRE  PARAMETERS 


Wire  length  in  meters 

0.1 

Wire  radius  in  meters 

0.0026 

Incident  angle  in  degrees 

90 

DG  pulse  minimum  10%  spread  in  nanoseconds 

0.28 

DG  pulse  maximum  10%  spread  in  nanoseconds 

0.30 

Time  window'  in  nanoseconds 

10.00 

Impact  delay  in  nanoseconds 

2 

Scattering  angle  in  degrees 

90 

Note  that  the  minimum  and  the  maximum  10%  spreads  determine  the  respective  a, 
and  a2  in  the  DG  waveform  of  equation  3.9. 

First,  the  frequency  spectrum  of  the  smoothed  impulse  response  was  calculated  from 
the  output  and  then  transformed  to  the  time  domain  by  an  FFT  [Ref.  18].  The  reason 
for  calculating  the  frequency  response  first  is  that  it  is  more  convenient  to  work  in  the 
frequency  domain  when  using  Raid's  method  to  extract  the  smoothed  impulse  response 
using  equation  3.6 

The  dotted  line  in  Figure  6  on  page  25  is  the  actual  measurement.  Because  the 
measurement  was  contaminated  by  noise,  the  result  did  not  agree  exactly.  The  contin¬ 
uous  line  represents  the  computed  thick  wire  smoothed  impulse  response  due  to  the 
same  excitation.  This  computation  was  performed  directly  in  the  time  domain  using  an 
E-field  integral  equation.  Measurements  for  other  types  of  structures  will  be  illustrated 
in  the  next  chapter. 


V.  EXPERIMENTAL  RESULTS 


A.  CANONICAL  MODEL  SCATTERING 

In  addition  to  the  10  cm  long  thick  wire  considered  in  Chapter  4,  four  additional 
canonical  targets  were  selected: 

•  10  centimeter  thin  wire  (broadside  and  canted  incidence  angles) 

•  8.1  centimeter  diameter  sphere 

•  12.2  centimeter  diameter  sphere 

•  15.5  centimeter  diameter  sphere 

The  Mie  series  program,  MIE,  was  used  to  compute  the  frequency  domain  scatter¬ 
ing  transfer  function  for  the  various  metallic  spheres.  These  transfer  functions  were  then 
used  as  inputs  to  the  SPRCT  program  to  compute  the  transient  scattering  responses  due 
to  the  specified  Double-Gaussian  (DG)  pulse  incident  plane  wave,  as  specified  by 
equation  3.9.  The  thin  wire  computations  were  performed  using  the  TDIG  program. 
This  program  utilizes  a  numerically  evolved  time-domain  electric  field  integral  equation 
for  the  induced  current  distribution  on  the  wire.  The  scattered  field  is  then  calculated 
by  appropriate  numerical  space-time  integrations  using  this  induced  current.  This  and 
the  sphere  scattering  programs  wrere  written  by  Morgan.  The  measured  and  computed 
DG  pulse  scattering  responses  are  overlayed  in  each  plot  to  validate  both  the  fidelity  of 
the  experimental  procedure  and  the  accuracy  of  the  deconvolution  algorithm.  The 
characteristics  of  each  canonical  target  scattering  case  are  summarized  in  Tables  9-13 
on  pages  28-33.  A  noteworthy  item  is  that  there  is  a  slight  bistatic  angle  (3  degrees) 
between  incident  and  scattered  aspects.  This  has  been  incorporated  into  the  computa¬ 
tions.  As  can  be  readily  observed,  the  agreement  between  measurements  and  computa¬ 
tions  in  Figures  7-11  on  pages  28-32  is  quite  good  in  all  cases.  The  additive  "rippling" 
of  the  experimental  results  is  due  in  part  to  the  residual  noise  and  clutter  contributions. 
This  latter  pollutant  includes  uncancelled  multipath  interactions  of  the  target  and 
canonical  sphere  scatterers  due  to  imperfect  absorption  by  the  anechoic  chamber  sur¬ 
faces.  Additional  differences  in  the  comparisons  are  due  to  the  approximations  which 
are  innate  to  the  computations. 


Figure  7.  The  12.2  cm  diameter  sphere  validation. 


Table  9.  12.2  CM  SPHERE  PARAMETERS 


Sphere  radius  in  meters 

0.061 

Distance  from  sphere  in  meters 

2.5 

Bistatic  angle  in  degrees 

3 

Scattering  plane 

E  plane 

Time  window  in  nanoseconds 

20.0 

Narrow  10%  pulse  width  in  nanoseconds 

0.28 

Wide  10%  pulse  width  in  nanoseconds 

0.30 

28 


Figure  8.  The  10  cm  thin  wire  validation  at  45  degrees 


Table  10.  45  DEGREE  THIN  WIRE  PARAMETERS 


Wire  length  in  meters 


Wire  radius  in  meters 


Incident  angle  in  degrees 


DG  pulse  minimum  10%  spread  in  nanoseconds 


DG  pulse  maximum  10%  spread  in  nanoseconds 


Scattering  angle  in  degrees 


No.  wire  segments 


0.00095 


45 


0.28 


0.30 


135 
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Figure  9.  The  10  cm  thin  wire  validation  at  90  degrees 


Table  1 1.  90  DEGREE  THIN  WIRE  PARAMETERS 


Wire  length  in  meters 

0.1 

Wire  radius  in  meters 

0.00095 

Incident  angle  in  degrees 

90 

DG  pulse  minimum  10%  spread  in  nanoseconds 

0.28 

DG  pulse  maximum  10%  spread  in  nanoseconds 

0.30 

No.  wire  segments 

20 

Scattering  angle  in  degrees 

90 

30 


MEASUREMENT  AND  COMPUTATION  COMPARISON 


8.1  CENTIMETER  DIAMETER  SPHERE 


Figure  10.  The  8.1  cm  diameter  sphere  validation 


Table  12.  8.1  CM  SPHERE  PARAMETERS 


Sphere  radius  in  meters 

0.0405 

Distance  from  sphere  in  meters 

2.5 

Bistatic  angle  in  degrees 

3.0 

Time  window  in  nanoseconds 

20.0 

Nyquist  frequency  in  Ghz 

12.775 

Narrow  10%  pulse  width  in  nanoseconds 

0.28 

Wide  10%  pulse  width  in  nanoseconds 

0.30 

MEASUREMENT  AND  COMPUTATION  COMPARISON 
1S.S  CENTIMETER  DIAMETER  SPHERE 


TIME  SERIES  IN  NANO  SECONDS 


Figure  11.  The  15.5  cm  diameter  sphere  validation 


B.  SCALE  MODEL  SCATTERING 

The  scattering  responses  of  the  wire  and  sphere  targets  due  to  the  incident  DG  in¬ 
cident  pulse  can  be  numerically  evaluated  using  straightforward  algorithms.  This  is  not 
the  case  for  more  complex  scattering  objects,  such  as  a  highly  realistic  model  of  an  air¬ 
craft.  Primitive  computational  algorithms  do  currently  exist  for  such  complex  structures 
but  the  poor  accuracy  of  their  results  would  be  insufficient  to  properly  validate  the 
measurements  being  made  here.  On  the  other  hand,  the  very  good  comparisons  shown 
in  the  previous  section  for  both  low-damped  (wires)  and  high-damped  (spheres)  strongly 


32 


i 


Table  13.  15.5  CM  SPHERE  PARAMETERS. 


Sphere  radius  in  meters 

0.0775 

Distance  from  sphere  in  meters 

2.5 

Bistatic  angle  in  degrees 

3.0 

Scattering  plane 

E  plane 

Time  window  in  nanoseconds 

20.0 

Nyquist  frequency  in  Ghz 

12.775 

No.  frequency  points 

257 

support  the  conjecture  that  high  quality  measurements  can  be  obtained  for  all  metallic 
targets  having  comparable  physical  dimensions. 

After  the  validations  were  made  using  wires  and  spheres,  extensive  measurements 
were  made  of  the  transient  backscattering  from  four  different  silver  coated  1/72  scale 
plastic  aircraft  models.  The  full-size  dimensions  of  each  aircraft,  whose  identities  will 
not  be  given  here,  are  shown  in  Tables  14-17  on  pages  34-37. 

Measurements  were  made  for  various  incident  aspects  in  the  wing- fuselage  plane  on 
each  model,  using  horizontal  incident  polarization:  0  degrees  (nose-on);  30  degrees;  90 
degrees  (broadside);  and  180  degrees  (tail-on).  In  addition,  a  broadside  look-down  as¬ 
pect  was  considered,  with  backscattering  measured  for  incident  linear  polarizations  both 
parallel  to  ("vertical  wing")  and  perpendicular  to  ("horizontal  wing")  the  fuselage. 
Samples  of  these  numerous  measurements  are  displayed  in  Figures  12-15  on  pages  34-37. 

The  complete  library  of  waveforms  is  being  used  to  support  research  in  demon¬ 
strating  the  accuracy  and  aspect  independence  of  natural  resonance  radar  target  iden¬ 
tification. 

Four  scale  model  aircraft  were  selected.  They  are  as  follows: 

•  Scale  model  1 

•  Scale  model  2 

•  Scale  model  3 

•  Scale  model  4. 

C.  SYSTEM  NOISE  ESTIMATION 

Noise  in  the  data  acquired  from  the  acquisition  program  was  investigated  in  this 
effort  in  order  to  develop  a  mechanism  to  estimate  the  signal-to-noise  ratio  (SNR).  The 
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Table  14.  DIMENSIONS  OF  SCALE  AIRCRAFT  MODEL  1 


Figure  13.  Scale  aircraft  model  1  at  000,  030,  090,  180  degree  aspect  angle 


Table  15.  DIMENSIONS  OF  SCALE  AIRCRAFT  MODEL  2 


Scale 

1/72 

Overall  length 

53  feet 

Extended  wingspan 

28  feet  7  inches 

Figure  14.  Scale  aircraft  model  comparisons  at  000  degree  aspect  angle 


Table  16.  DIMENSIONS  OF  SCALE  AIRCRAFT  MODEL  3 


Scale 

1/72 

Length 

56  feet 

Wingspan 

40  feet  8  inches 

Height 

15  feet  4  inches 

Figure  15.  Scale  aircraft  model  comparisons  090  degree  aspect  angle 


Table  17.  DIMENSIONS  OF  SCALE  AIRCRAFT  MODEL  4 


Scale 

1/72 

Length 

63.8  feet 

Wingspan 

42.8  feet 

Height 

18.6  feet 

37 


SNR  is  a  major  determining  factor  for  the  measurement  fidelity.  Much  of  this  SNR 
assessment  involved  the  computing  some  statistics  of  the  waveform.  All  of  these  statis¬ 
tical  calculations  were  performed  using  GRAFSTAT  with  APL. 

In  the  acquisition  step,  the  calibration,  the  background,  the  target  and  a  second 
background  were  measured  consecutively.  The  two  backgrounds  were  measured  for  the 
purpose  of  noise  estimation.  In  the  following  analysis,  a  2048  ensemble-average  mode 
was  used  to  acquire  the  data.  The  acquisition  of  each  waveform  took  about  9  minutes 
to  complete. 

In  a  separate  investigation,  the  stationarity  of  the  system  noise  was  investigated  by 
acquiring  seven  consecutive  background  measurements.  These  backgrounds  were  then 
subtracted  from  each  other  to  observe  the  noise  occuring  during  the  consecutive  meas¬ 
urement.  For  example,  the  second  background  and  third  background  were  subtracted, 
as  is  shown  in  Figure  16  on  page  39  along  with  various  statistical  parameters  in  Table 
18  on  page  39  mentioned,  the  statistics  of  each  subtracted  waveform  was  calculated  by 
the  GRAFSTAT  software  packgagc.  The  variance  of  the  noise  can  be  calculated  from 
the  standard  deviation  and  expected  power  can  be  approximated  by  this  variance  when 
the  noise  has  a  small  mean.  This  approximation  was  reasonable  because  the  squares  of 
the  mean  for  the  actual  measurement  data  are  very  small,  (  sec  Figure  16  on  page  39  ). 
The  computed  noise  power  is  plotted  versus  time  in  Figure  17  on  page  40.  As  can  be 
seen,  the  measured  noise  power  varied  only  slightly  with  time.  A  linear  least-squares  lit 
yields  the  equation 

Y  »  -93.042  +  0.035727 A'  ,  (5.1) 

where  Y  is  the  expected  noise  power  and  X  is  the  time  in  minutes.  This  model  was 
based  on  only  seven  background  measurements  with  only  2048  averages  per  measure¬ 
ment.  For  practical  purposes  the  system  noise  can  be  assumed  to  be  fixed  in  time,  de¬ 
pending  only  upon  the  number  of  ensemble  averages  used  to  complete  the  waveform 
measurement. 

To  calculate  the  SNR,  the  thick  wire  measurement  was  used.  The  signal  plus  noise 
waveform  is  displayed  in  Figure  18  on  page  41  while  the  computed  statistical  parameters 
are  given  in  Table  20  on  page  41. 

This  waveform  was  time  windowed  to  eliminate  the  noise  outside  of  the  window' 
before  the  statistical  parameters  were  calculated.  The  signal  plus  noise  power  was  cal¬ 
culated  by  way  of  the  variance  using  the  following  equation. 
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EXAMPLE  NOISE  WAVEFORM 


Figure  16.  Noise  waveform  statistics 


Table  18.  STATISTICS  OF  NOISE  WAVEFORM 


Mean 

-3.021  IE-10 

Std.  deviation 

0.000023345 

Median 

-7.019E-8 

Maximum 

0.000066731 
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NOISE  ANALYSIS 


NMK  ran  OP  W  TAIL 


Figure  17.  Noise  power 


SCATTER  PLOT  TABLE 


X 

TIME 

Y 

10x10*P0WERS 

SELECTION 

ALL 

X  LABEL 

TIWC  IN  MINUTE 

Y  LABEL 

08 

NO.  OF  ELEKNTS 

S 

CORRELATION  XY 

0.90183 

RK  CORRELATION 

0.93779  T-1.696 

X  UEAN 

31.9 

STD.  DEVIATION 

19.837 

9-PERCENTILE 

9 

25-PERCENTILE 

18 

MEDIAN 

31.9 

75-PERCENTILE 

45 

95-PERCENTILE 

94 

X  MIN. 

9  18  27 

X  MAX. 

94  45  36 

Y  MEAN 

-91.917 

STD.  DEVIATION 

0.99986 

9-PERCENTILE 

-92.79 

29-PERCENTILE 

-92.716 

kCDIAN 

-92.343 

79-PERCENTILE 

-90.889 

95-PERCENTILE 

-90.869 

Y  MIN 

-92.78  -92.718  -92.838 

Y  MAX 

-90.889  -90.689  -92.049 

Table  19.  STATISTICS  OF  NOISE  MODEL 


Mean 

-91.917 

Std.  deviation 

0.99986 

Median 

-92.343 

Maximum 

-90.669 

THICK  WIRE  SIGNAL  PLUS  NOISE 
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BATTER  PUTT  TABLE 
X  :tXt 

V  :TOTAL*<TTMHX-rmrBT) 

(ELECTION  :ALL 

X  LABEL  :TUC  IN  NANO 

V  LABEL  sAtPUTVJOC  OA  HAND 


NO.  OF  ELEMENT! 

>312 

CORRELATION  XV 

:0.00BMB1 

IK  CORRELATION 

>0.00*1*07  T-C.1JN1 

X  MEAN 

!10 

ITS.  DEVIATION 

>*.7104 

B-FERCENT1LE 

>0.*823* 

2S-PERCEXTILE 

K.noz 

MEDIAN 

1 10 

76-PBCEXTILE 

>11.01 

•8-PEHCENTILE 

>11.01! 

X  MIN. 

:0  0. 03*131  0.07*27* 

X  MAX. 

>20  1t.N1  19.822 

V  ICAN 

:**.2N1E-* 

ITS.  DEVIATION 

>0.00033*2 

6-PEJCENTILE 

>-0.0000*722 

26-PERCENTILE 

>0 

MEDIAN 

>0 

7S-PERCEMTILE 

>0 

■S^ERCENTILE 

>0.0006*718 

V  MIN 

>-0. 0014*22  -0.00130* 

V  MAX 

>0.001412*  0.00131*  0. 

Figure  18.  Signal  plus  noise  waveform  statistics 


Table  20.  STATISTICS  OF  THICK  WIRE  SIGNAL  PLUS  NOISE 


WAVEFORM 

Mean 

-6.299  IE-6 

Std.  deviation 

0.0003362 

Median 

0 

Maximum 

0.0014128 

(5.2) 


E(X2)  =  V ariance  +  E2(X) 

where  E(X)  is  mean  of  the  noise  and  E{Xi)  is  the  mean  square  value  of  the  noise  (noise 
power).  Note  that  the  second  term  on  the  right  side  of  equation  5.2  was  close  to  zero, 
so  it  can  be  neglected.  As  a  result,  the  expected  power  is  equal  to  the  variance.  The 
signal  plus  noise  power  was  equal  to  0.0003382*  or  -69.41  dB.  The  noise  expected  power 
was  calculated  by  equation  5.1  and  is  given  by  -92.64  dB  for  the  average  time  of  10 
minutes.  Since  the  noise  power  was  so  small  compared  to  the  signal  plus  noise  power, 
the  subtraction  between  the  signal  plus  noise  and  the  noise  was  not  necessary.  The 
noise  power  when  time  windowed  was  reduced  by  a  factor  of  200/512,  or  4.08  dB.  The 
result  for  the  SNR  is  27.31  dB  in  the  case  of  the  thick  wire  measurement.  Note  that  this 
SNR  is  about  10  dB  better  than  that  provided  for  the  same  target  using  the  Tektronix 
equipment,  as  carried  out  by  McDaniel  (Ref.  1 1). 


VI.  CONCLUSIONS 


A.  SUMMARY 

This  thesis  describes  the  theory  and  implementation  of  a  newly  updated  transient 
scattering  measurement  facility.  An  overview  of  the  historical  precedents  and  moti¬ 
vations  for  this  development  was  given  in  Chapter  I.  This  was  followed  in  Chapter  II 
by  a  detailed  description  of  the  new  equipment  being  used  as  well  as  the  physical  struc¬ 
ture  of  the  anechoic  chamber. 

Building  upon  this  overview  of  the  hardware  aspects,  the  conceptual  and  math¬ 
ematical  modeling  of  the  scattering  range  is  undertaken  in  Chapter  III.  Using  the 
physical  model  of  the  radiating  and  receiving  antennas  in  the  anechoic  chamber,  the 
various  electromagnetic  interactions  between  the  scattering  target  and  the  chamber  are 
represented  by  a  comprehensive  linear  system  model.  The  topology  of  the  various  fre¬ 
quency  domain  transfer  functions  in  this  model  is  then  used  to  justify  the  measurement 
and  pre-processing  steps  required  to  estimate  the  smoothed  scattering  impulse  response, 
due  to  a  specified  double-Gaussian  incident  plane  wave.  There  are  three  measurements 
required  for  each  target  impulse  response  estimation.  Aside  from  the  actual  target 
measurement,  a  background  measurement  is  needed  to  subtract  out  the  directly  received 
signal  from  the  transmitting  antenna,  in  addition  to  that  scattered  from  the  various 
chamber  surfaces.  A  third  measurement,  for  a  metallic  sphere,  is  needed  to  allow  for 
compensation  of  the  transfer  function  characteristics  of  the  transmitting  and  receiving 
antennas  in  the  chamber.  By  effecting  a  deconvolution  comparison  with  a  highly  accu¬ 
rate  numerical  calculation  of  the  transient  scattering  response  of  the  sphere  target,  this 
final  measurement  also  permits  the  estimation  of  the  target's  response  to  the  same 
excitation  as  is  used  in  the  sphere  computation. 

Chapter  IV  considers  the  details  of  the  actual  software  that  performs  the  signal  ac¬ 
quisition,  signal  averaging  and  deconvolution  post-processing  of  the  measured 
waveforms.  The  acquisition  and  signal  averaging  software  is  a  blend  of  a  commercial 
IEEE-488  bus  controller  (GURL-II)  and  a  custom  driver,  designed  as  part  of  this  thesis 
effort.  The  deconvolution  algorithm  provides,  first  of  all,  graphically  displayed  iterative 
background  subtractions  from  both  the  target  and  sphere  measurements  to  allow  com¬ 
pensation  for  time  drifts.  Secondly,  the  final  deconvolution  is  performed  via  a  variant 
on  Riad's  method,  as  discussed  in  Chapter  III. 


The  validations  shown  in  Chapter  V  speak  for  themselves  as  to  the  expected  meas¬ 
urement  quality  obtained  using  this  new  transient  scattering  range.  Both  low-Q  and 
high-Q  type  targets  are  considered  by  the  respective  comparisons  of  sphere  and  thin-wire 
scattering  measurements  with  computed  smoothed  impulse  responses.  Further  illus¬ 
trations  of  the  type  of  measurements  to  be  performed  with  this  range  are  shown  in  the 
scale  model  aircraft  backscattering  signatures.  A  simple  software  modification  will  allow 
the  display  of  very  broadband  frequency  domain  radar  cross  sections  (RCS)  for  scale 
model  targets. 

To  quantitatively  investigate  the  fidelity  of  the  measurement  system,  a  noise  esti¬ 
mation  of  the  new  implementation  was  investigated  and  was  modelled.  The  particular 
signal-to-noise  ratio  (SNR)  of  a  thick  wire  scattering  measurement  was  calculated.  A 
27  dB  SNR  of  this  case  demonstrated  the  relatively  high  fidelity  of  the  measurement 
system  for  the  case  small  RCS  targets,  including  silver  coated  1/72  scale  model  aircraft. 

B.  FUTURE  CONSIDERATIONS 

Even  though  the  demonstrated  system  performance  was  excellent,  there  is  no  such 
thing  as  too  high  an  SNR  in  scattering  measurements.  In  this  regard,  the  power- 
bandwidth  product  has  the  potential  to  be  enhanced  by  future  efforts. 

The  1  to  7  GHz  passband  of  the  scattering  system  was  limited  by  the  transmitting 
GaAS  power  amplifier.  On  the  other  hand,  the  bandwidth  of  the  HP  DPO  system  can 
go  as  high  as  20  Ghz,  while  the  current  antennas  have  a  1  to  12  GHz  passband.  Future 
use  of  broader  bandwidth  antennas  and  power  amplifier,  having  a  rated  RMS  output 
power  of  2  watts  or  more,  would  increase  the  overall  system  performance.  More  accu¬ 
rate  smoothed  impulse  responses  could  be  achieved,  as  well  as  wider  bandwidth  RCS 
measurements. 

Because  of  the  limited  scope  of  this  thesis,  the  noise  power  estimation  was  only  an 
approximation.  A  more  comprehensive  noise  analysis  should  be  carried  out  in  order  to 
quantify  and  minimize  the  noise  power  sources. 

The  final  consideration  is  the  future  investigation  of  alternate  mechanisms  for  pro¬ 
viding  broadband  scattering  measurements.  One  such  mature  technique  is  through  a 
stepped  frequency  continuous  wave  system,  employing  a  network  analyzer  front-end  as 
the  coherent  receiver,  perhaps  preceded  by  a  low-noise  amplifier  to  reduce  the  overall 
system  noise  figure.  A  somewhat  radical  technique,  one  that  is  currently  being  investi¬ 
gated  at  N'PS,  is  to  use  a  very  broadband  amplified  noise  source  as  the  transmitting 
source,  while  employing  the  HP  DPO  as  a  two-channel  correlation  receiver.  This 
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method  is  akin  to  a  dynamic  matched  filter,  offering  enhanced  bandwidth,  at  relatively 
low  cost  vis-a-vis  the  stepped  frequency  setup.  Additionally,  the  potential  exists  for  both 
tactical  and  strategic  employment  of  such  random  source  systems,  either  as  active 
"pseudo-jammers"  or  as  totally  passive  radars  relying  upon  either  natural  or  man-made 
random  signal  sources  for  the  illumination  of  targets. 
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APPENDIX  A.  THE  IMPULSE  RESPONSE  DECONVOLUTION 


The  following  computer  program  calculates  the  Impulse  Response  from  the  meas¬ 
urement  data  as  described  in  chapter  4.  The  program  is  written  in  RM, FORTRAN. 

The  subroutine  PL08  and  Pll  require  the  GRAFMATIC  routine  library.  The  pro¬ 
cedure  to  compile,  link  and  run  this  source  code  refers  to  the  RM/FORTRAN  User's 
Guide  Version  2.1!  (DOS). 


C  THIS  IS  THE  IMPULSE  RESPONSE  DECONVOLUTION  * 
C  WRITTEN  BY  LT.  JG  SOONPUEN  SOMAPEE  Feb  88  * 
C  ADAPTATION  OF  ORIGINAL  TEKTRONIX  PROGRAM  * 
C  BY  PROF.  M.  A. MORGAN  Circa  1983-85  * 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


WHERE*  * 

DATE  IS  THE  DATE  OF  MEASUREMENT  * 

DATAB  IS  THE  NUMBER  OF  DATA  BLOCKS  * 

FI  IS  THE  DUMMY  STRING  * 

ID  IS  THE  IDENTIFICATION  STRING  * 

K  IS  THE  NUMBER  OF  POINTS  * 

L9  IS  THE  RELATIVE  SMOOTHING  * 

MAG  IS  THE  MAGNITUDE  OF  Xo*Ht(f)  * 

NAME4  IS  THE  NAME  OF  THE  TARGET  DATA  * 

N  IS  THE  NUMBER  OF  POINTS  * 

NUMPT  IS  THE  SAME  AS  K  * 

NEXPT  IS  INDICATEES  POWER-OF-2  EXPONENT  * 
WILL  BE  GENERATE  BY  PREFFT  * 

SET  TO  -1  INDICATES  ERROR  * 

R  IS  THE  REAL  PART  OF  THE  INVERSE  * 

FOURIER  TRANSFORM  OF  Xo*Ht(f)  * 

SUBAV  IS  THE  NUMBER  OF  SUB-AVERAGES  * 

T  IS  THE  TIME  INTERVAL  * 

TIME  IS  THE  TIME  INTEVAL  * 

TITLE  IS  THE  TITLE  NAME  FOR  PLOTTING  * 

TUT  IS  THE  FREQUENCY  DOMAIN  OF  Xo*H(f)* 

TYPE  IS  THE  WAVEFORM  TYPE  * 

W  IS  THE  COMPLEX  EXPONENTIAL  ARRAY  * 

GENERATED  BY  PREFFT  SUBROUTINE  * 

Xo*Ht( f)  IS  THE  GAUSSIAN  WEIGHTED  COSINE  * 

XS  IS  THE  STRING  FOR  PLOTTING  ROUTINE* 

FOR  THE  HEADER  * 

DRIVER  * 

Y4  IS  TT  CALIBRATION  SPHERE  * 

SUBTRACTED  WAVEFORM  (FREQUENCY)  * 

* 

Y5  IS  THE  TARGET  SUBTRACTED  WAVEFORM  * 

FREQUENCY  DOMAIN  * 

* 


C  YSC  IS  THE  COMPUTED  SCATTERED  ELECTRIC* 
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nonooono 


Y1 

FIELD  FOR  THE  CANONICAL  TARGET 

IS  THE  CALIBRATION  MEASUREMENT 

•k 

* 

Y2 

IS 

THE 

BACKGROUND  MEASUREMENT 

* 

Y3 

IS 

THE 

TARGET  MEASUREMENT 

* 

YS 

IS 

THE 

SAME  AS  XS 

ie 

YN 

IS 

THE 

YES  OR  NO  STRING 

* 

XMI1 

IS 

THE 

MINIMUM  VALUE 

* 

XMA1 

IS 

THE 

MAXIMUM  VALUE 

* 

CHARACTER* 64  TITLE , XS , ID , YS , DATE , TYPE , NAME 
CHARACTER*64  TYPE 1 , TYPE  2 , TYPE  3 , F 1 
CHARACTER* 16  NAME 1 .NAME 2 .NAME 3 
CHARACTER*64  ID1.ID2.ID3 
CHARACTER* 1  YNS.YNS1 

COMPLEX  YSC( 1025) ,Y4( 1025) ,Y5( 1025) ,W( 1025) 
REAL  T,T2  .TIME 

REAL  Yl( 1024), Y2( 1024) ,Y3( 1024) 

INTEGER  N , SUBAV , DATAB , NUMPT 


YNS='N' 


READ  DATA  FOR  COMPUTED  SPHERE  THEN 
CALLING  SUBROUTINE  FFF  TO  DRAW  THE  TIME 
SERIES  AND  SPECTRAL  PLOT 
■  *** . 


WRITE(*,*) 

WRITE(*,*) 


WRITE(*,*) 

'* 

DECONVOLUTION  PROGRAM 

BY  DR.  M.  A.  MORGAN 

WRITE(*,*) 

'* 

ADAPTED 

TO  IBM  BY  LT.  JG. 

SOONPUEN  SOMAPEE 

*' 

WRITE(*,*) 

WRITE(*,*) 

'* 

FORMAT  OF  DATA  FILE: 

WRITE(*,*) 

WRITE(*,*) 

»* 

1. 

READ 

ID 

(FORMAT  CHARACTER*64 ) 

*' 

WRITE(*,*) 

'* 

2. 

READ 

TYPE 

(FORMAT  CHARACTER*64) 

WRITE(*,*) 

3. 

READ 

DATE 

(FORMAT  CHARACTER*64 ) 

*' 

WRITE(*,*) 

'* 

4. 

READ 

FI 

(FORMAT  CHARACTER*64 ) 

WRITE(*,*) 

'* 

5. 

READ 

SUB -AVERAGE 

(FORMAT  INTEGER) 

WRITE(* ,*) 

6. 

READ 

TIME  WINDOW 

(FORMAT  REAL(F5.  2)) 

*• 

WRITE(*,*) 

7. 

READ 

DATA  BLOCK 

(FORMAT  INTEGER) 

*' 

WRITE(*,*) 

'* 

8. 

READ 

NUMBER  OF  POINTS 

(FORMAT  INTEGER) 

WRITE(*!*)  'FILE  NAME  FOR  COMPUTED  CALIBRATION  SPHERE: 
READ(*, 173)  NAME 


YSC  =  COMPUTED  SPHERE  COMPLEX  FREQUENCY  DOMAIN  DATA 
NUMPT  =  NUMBER  OF  POINTS  IN  THE  OUTPUT  OF  THIS  ROUTINE 
TIME  =  OUTPUT  TIME  INTERVAL 


CALL  FFFKNAME, YSC, W.NEXPT, TIME) 


9110  CONTINUE 


C  READ  THE  DATA  FILES  FOR  CALIBRATE  AND  BACKGROUND 
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THEN  OVERLAY  THESE  ONTO  THE  SAME  AXIS. 

PLOT77  ALLOWS  THE  OPERATOR  TO  SHIFT  ONE  WAVEFORM 
PRIOR  TO  SUBTRACTING  THE  BACKGROUND. 


C 
C 
C 
C 

C  THE  RESULT  FROM  SUBTRACTION  IS  TRANSFORMED  TO 
C  FREQUENCY  DOMAIN  AND  SAVED  TO  THE  COMPLEX  ARRAY  Y4 
C 

C  NAME1  =  NAME  OF  THE  FILE  TO  SENT  TO  PLOT 7 7  ROUTINE 
C  (CALIBRATE  DATA) 

C  NAME 2  =  NAME  OF  THE  BACKGROUND  DATA 
C  ID1  =  IDENTIFICATION  CHR$  FOR  CALIBRATE 

C  Y4  =  COMPLEX  ARRAY  OUTPUT  (CALIBRATE  ) 

C***  *********  A1  ft  ************  irkirkirkirirklrlrirk  ■*'*  ***********  A1**1**** 

WRITE(*,*)  1  **********  * A A *******  A************ A******* 1 
WRITE(*,*)  'FILE  NAME  FOR  CALIBRATION  SPHERE:  ’ 
READ(*,173)  NAME1 

CALL  INPP(NAME1 , ID1 ,TYPE1 .DATE ,F1 , SUBAV.T2.DATAB .NUMPT.Yl) 
WRITE(*,*)  '  **************  *  *  ******  **  ***  A************* 1 
WRITE(*|*)  'FILE  NAME  FOR  SPHERE  BACKGROUND  WAVEFORM:  ' 
READ(*,  173)  NAME 2 


CALL  INPP( NAME 2 , ID2 , TYPE  2 , DATE , F 1 , SUBAV , T2 , DATAB , NUMPT , Y2 ) 
CALL  PLC77(Y1,Y2,T2, NUMPT, ID1) 


9112  CONTINUE 

WRITE(*,*)  '  ***************************************** 1 
WRITE(*i*)  'FILE  NAME  FOR  TARGET  WAVEFORM:  ' 

READ(*, 173)  NAME 3 

C*********************  *****  A  *  ***  *  *  *  ************** 


CALL  I NPP ( NAME 3.ID3, TYPE 3 , DATE , F 1 , SUB AV , T2 , DATAB , NUMPT , Y3 ) 


WRITE(*,*) 

WRITE(*,*)  'USE  PREVIOUS  BACKGROUND  FOR  THE  TARGET  7  (Y/N): 
READ(*, 173)  YNS1 

IF(YNS1.  EQ.  jY'.OR.YNSl.EQ.  'y')  GO  TO  9222 

WRITE(*!*)  'FILE  NAME  FOR  TARGET  BACKGROUND  WAVEFORM:  ’ 
READ(*, 173)  NAME 2 


I 


CALL  INPP ( NAME2 , ID2 , TYPE  2 , DATE , F 1 , SUB AV , T2 , DATAB , NUMPT , Y2 ) 

9222  CONTINUE 
C 

C  REPEATING  SHIFT  AND  SUBTRACT  OPERATIONS  ON  TARGET  AND  BACKGROUND 
C  NAME1=  TARGET  DATA  FILE 

C  NAME 2=  BACKGROUND  DATA  FILE 

C  OUTPUT 

C  ID1  =  IDENTIFICATION  OF  THE  TARGET 

C  Y5  =  COMPLEX  ARRAY  OF  FREQUENCY  DOMAIN  TARGET  DATA 

Q******************************* *  A  **  *  *****  A******* 


CALL  PL077(Y3,Y2,T2, NUMPT, ID3) 

DO  987  1=1, NUMPT 

Y4( I)=CMPLX(Y1( I) ,0.  0) 

Y5( I)=CMPLX(Y3( I) ,0.  0) 

987  CONTINUE 

MODE=0 

CALL  FFT( NUMPT, MODE, TIME, NEXPT,W,Y4) 


CALL  FFTC NUMPT, MODE, TIME, NEXPT,W,Y5) . 

CALL  DECON3C NUMPT, YSC , Y4 , Y5 , I D3 , TYPE 3 , DATE , F 1 , SUB AV , T2 , DATAB ) 

WRITEC*’*)  'WOULD  YOU  LIKE  TO  PROCESS  ANOTHER  TARGET  ?  (Y/N):  ' 
READ(*, 173)  YNS 

IF(YNS.  EQ.  '  Y'  .OR.  YNS.  EQ.  V  )  GO  TO  9112 
WRITEC*,*)  1  ********************** k a ****★*★★* 1 
WRITEC*,*)  '*  eeeee  eeeeee  e  *' 

WRITEC*,*)  '*  eee  eee  eeee  *' 

WRITEC*,*)  ' eeeeeeeeeeeeeeeeeeeeeeeeeeeee  *' 

WRITEC*,*)  '  eeeeeeeeeeeeeeeeeeeeeeee  *' 

WRITEC*,*)  '  eeeeeeeeeeeeeeeeeeeeee  *' 

WRITEC*,*)  1  ****^^*******<r»*****^********* 1 
WRITEC*,*)  ' ********************************* * 

173  FORMAT(A) 

STOP 

END 


SUBROUTINE  INPP(  NAME, ID, TYPE, DATE, FI, SUBAV.T2, DATAB, NUMPT.Y) 


C 

C  INPUT  NAME 

C  OUTPUT  ID, TYPE, DATE, FI, SUBAV.T2, DATAB, NUMPT.Y 

C  NAME  =  CHR$  NAME  OF  THE  DATA  FILE 

C  ID  =  CHR$  TARGET  IDENTIFICATION 

C  TYPE  =  CHR$  WAVE  FORM  IDENTIFICATION 

C  DATE  =  CHR$  DATE 

C  INTEGER 

C  SUBAV  =  NUMBER  OF  SUB -AVERAGES  USED  TO  ACQUIRE  THE  DATA 
C  DATAB  =  DATA  BLOCK 

C  FI  DUMMY 

C  NUMPT  =  NUMBER  OF  POINTS 

C  REAL 

C  T2  TIME  WINDOW 

C  Y  =  Y  ARRAY  OUTPUT 


CHARACTER* 16  NAME 
CHARACTER*64  ID, FI 
CHARACTER*64  TYPE, DATE 
REAL*4  Y( 1500) 

INTEGER  SUBAV, DATAB, NUMPT, K 

PI=3.  141592654 

K=1 

OPEN  (UNIT=10 ,FILE=NAME , STATUS3' OLD' ) 
C  READ  THE  HEADER 

READ( 10,5)  ID 
READ( 10,5)  TYPE 
READ( 10,5)  DATE 
READ( 10,5)  FI 
READC 10,2030)  SUBAV 
READ( 10,2010)  T2 
READC 10, 2030)  DATAB 
READ( 10,2030)  NUMPT 
10  READC 10, 2000, END=20)  YCK) 
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K=K+1 


17 

GOTO  10 

20 

CLOSE(IO) 

2000 

FORMAT(E12.  6) 

2010 

FORMATC F5.  2) 

2030 

F0RMAT(I5) 

5 

FORMATC A) 

K=x- 1 

IF(K.  EQ.  NUMPT)  GOTO  3000 

WRITE(*,*)  'CHECK  THE  DATA  -  if  POINTS  MISMATCH* 
STOP 

3000  RETURN 
END 


SUBROUTINE  FFF1(  NAME, YY,W,NEXPT, TIME) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


WHERE  * 

DATE  IS  THE  DATE  OF  MEASUREMENT  * 

DATAB  IS  THE  NUMBER  OF  DATA  BLOCKS  * 

FI  IS  THE  DUMMY  STRING  * 

ID  IS  THE  IDENTIFICATION  STRING  * 

K  IS  THE  NUMBER  OF  POINT  * 

L9  IS  THE  RELATIVE  SMOOTHING  * 

MODE  IS  EQUAL  TO  0  MEAN  FORWARD  * 

FOURIER  TRANSFORM 

MAG  IS  THE  MAGNITUDE  OF  Xo*Ht(f)  * 

NAME4  IS  THE  NAME  OF  THE  TARGET  DATA  * 

N  IS  THE  NUMBER  OF  POINTS  * 

NUMPT  IS  THE  SAME  AS  K  * 

NEXPT  IS  THE  POWER-OF-2  EXPONENT  * 

WILL  BE  GENERATED  BY  PREFFT  * 

SET  TO  -1  INDICATES  ERROR  * 

R  IS  THE  REAL  PART  OF  THE  INVERSE  * 

FOURIER  TRANSFORM  OF  Xo*Ht(f)  * 

SUBAV  IS  THE  NUMBER  OF  SUB  AVERAGE  * 

T  IS  THE  TIME  INTERVAL  * 

TIME  IS  THE  TIME  INTEVAL  * 

TITLE  IS  THE  TITLE  NAME  FOR  PLOTTING  * 

TUT  IS  THE  FREQUENCY  DOMAIN  OF  Xo*H(f)* 

TYPE  IS  THE  WAVEFORM  TYPE  * 

W  IS  THE  COMPLEX  EXPONENTIAL  ARRAY  * 

GENERATE  BY  PREFFT  SUBROUTINE  * 

Xo*Ht( f )  IS  THE  GAUSSIAN  WEIGHTED  COSINE  * 

XS  IS  THE  STRING  FOR  PLOTTING  ROUTINE* 

FOR  THE  HEADER  * 

DRIVER  * 

Y4  IS  THE  CALIBRATION  SPHERE  * 

SUBTRACTED  WAVEFORM  (FREQUENCY)  * 

* 

Y5  IS  THE  TARGET  SUBTRACTED  WAVEFORM  * 

FREQUENCY  DOMAIN  * 

* 

YSC  IS  THE  COMPUTED  SCATTERED  ELECTRIC* 

FIELD  FOR  THE  CANONICAL  TARGET  * 

Y1  IS  THE  CALIBRATION  MEASUREMENT  * 

Y2  IS  THE  BACKGROUND  MEASUREMENT  * 
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oooooooooo  o a n o o n o 


Y3 

YS 

YN 


XHI1 

XMA1 


IS  THE  TARGET  MEASURED  * 
IS  THE  SAME  AS  XS  * 
IS  THE  YES  OR  NO  STRING  * 
METER  i 
IS  THE  MINIMUM  VALUE  * 
IS  THE  MAXIMUM  VALUE  * 


:**********  A  A**  A***  ■kirlrkirititir/rk'k'kirA  A  A*'*'***  *******  "*'*"<. 


CHARACTER* 16  NAME 

CHARACTER*64  ID , TYPE , DATE , TITLE , XS , YS , F 1 
COMPLEX  W( 1025 ),YY( 1025) 

REAL  T 

REAL  Y( 1024) ,TIME ,T2 ,XMINS ,MAG( 1024) ,XMAG 
INTEGER  NUMPT , DATAB , SUBAV , MODE , N 

C*+***-k*1'-kick****it+****-ir****-k*'kir>'*irk*******1'**k1'*i' 

CALL  INPP(  NAME, ID, TYPE, DATE, FI, SUBAV.T2, DATAB, NUMPT, Y) 
TIME=T2 / ( FLOATC  NUMPT- 1 )) * IE - 9 
XMINS=0. 0 
N=NUMPT 

CALL  PRETIM( TITLE, XS.YS) 

CALL  P 1 1( TITLE ,N,XMINS,T2,Y,XS,ID,YS) 

C  SET  MODE  EQUAi  TO  0  TO  FIND  FOURIER  SERIES 
MODE=0 

CALL  PREFFT  (N,MODE,NEXPT,W) 

DO  1234  I  =1  , NUMPT 

YY(I)=CMPLX(  Y( I) ,0.  0  ) 

1234  CONTINUE 

T=TIME 

CALL  FFT  (N,MODE ,T,NEXPT,W,YY) 

C  CHANGING  THE  OUTPUT  YY  TO  MAGNITUDE  AND  PHASE 
N=NUMPT 

CALL  PREFRE( TITLE, XS,YS) 

CALL  MAGNITUD(YY,N,MAG) 

XMAG=1.  O/TIME/2.  0*lE-9 

CALL  P11(TITLE ,N,XMINS ,XMAG,MAG,XS ,ID,YS) 

RETURN 
END 


* 


FFT  PROGRAM  FROM  THE  DIGITAL  SPECTRAL  ANALYSIS  * 
WITH  APLICATIONS  PRENTICE-HALL, INC. .  * 

ENGLEWOOD  CLIFF, NEW  JERSEY  07632  1987  * 

S.  LAWRENCE  MARPLE  * 

* 


*  jHHHt  A  A1*'****  A  A1  AA'****AAr**,*'A"Ar**iA**lA  A  A  AA'A'ArifrArArArArA'A'A'A  A  A  A'AAArATAA'ArAr 

SUBROUTINE  PREFFT  (N.MODE ,NEXP,W) 


Input  Parameters: 

N  -  Number  of  data  samples  to  be  processed  (integer-must  be 
power  of  two) 

MODE  -  Set  to  0  for  discrete-time  Fourier  series  (Eq.  2.C. 1)  o 
1  for  inverse  (Eq.  2.C.  2) 


Output  Parameters: 


oononooooonononno  o  10  t_.  m  cn  000000000 


I 


NEXP  -  Indicates  power-of-2  exponent  such  that  N=2**NEXP  . 

Will  be  set  to  -1  to  indicate  error  condition  if  N 
is  not  a  power  of  2  (this  integer  used  by  sub.  FFT) 

W  -  Complex  exponential  array 

Notes: 

External  array  W  must  be  dimensioned  .GE.  N  by  calling  program. 

COMPLEX  W( 1024) ,C1,C2 

NEXP=1 

NT=2**NEXP 

IF  (NT  .GE.  N)  GO  TO  10 

NEXP=NEXP+1 

GO  TO  5 

3  IF  (NT  .EQ.  N)  GO  TO  IS 
NEXP=-1 
RETURN 

5  S=8. *ATAN( 1.  )/FL0AT(NT) 

C1=CMPLX(C0S( S) , -SIN( S) ) 

IF  (MODE  .NE.  0)  C1=C0NJG(C1) 

C2=(  1.  ,0.) 

DO  20  K=1,NT 
W(K)=C2 

0  C2=C2*C1 

RETURN 
END 

*******************************  A  *  ****** **** ***** 


SUBROUTINE  FFT  (N, MODE, T, NEXP, W,X) 

Input  Parameters: 

N, MODE, NEXP, W  -  See  parameter  list  for  subroutine  PREFFT 
T  -  Sample  interval  in  seconds 

X  -  Array  of  N  complex  data  samples,  X(l)  to  X(N) 

Output  Parameters: 

X  -  N  complex  transform  values  replace  original  data  samples 
indexed  from  k=l  to  k=N,  representing  the  frequencies 
(k-l)/NT  hertz 

Notes: 

External  array  X  must  be  dimensioned  .GE.  N  by  calling  program. 

COMPLEX  X( 1024) ,W( 1024), Cl, C2 

MM=1 

LL=N 

DO  70  K=1 ,NEXP 
NN=LL/2 
JJ=MM+1 

DO  40  1=1 ,N,LL 
KK=I+NN 
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C1=X(I)+X(KK) 

X(KK)=X( I)-X(KK) 

40  X(I)=C1 

IF  (NN  .EQ.  1)  GO  TO  70 
DO  60  J=2,NN 
C2=W( JJ) 

DO  50  I=J,N,LL 
KK=I+NN 
C1=X(I)+X(KK) 
X(KK)=(X( I) -X(KK) )*C2 
50  X(I)=C1 


60 

JJ=JJ+MM 

LL=NN 

MM=MM*2 

70 

CONTINUE 

NV2=N/2 

NM1=N-1 

J=1 

DO  90  1=1, NM1 

IF  (I  .GE.  J) 
C1=X( J) 
X(J)=X(I) 
X(I)=C1 

GO  TO  80 

80 

K=NV2 

85 

IF  (K  .GE.  J) 
J=J-K 

K=K/2 

GO  TO  85 

GO  TO  90 

90 

J=J+K 

IF  (MODE  .EQ.  0) 

S=T 

IF  (MODE  .NE.  0) 
DO  100  1=1, N 

S=l.  /(T*FLOAT(N)) 

100 

X(I)=X(I)*S 

RETURN 

END 


SUBROUTINE  PRETIM( TITLE, XS,YS) 
CHARACTER*64  TITLE, XS.YS 


TITLE=’TIME  SERIES  PLOT' 

XS  =' X-AXIS  TIME  IN  NS  ' 

YS  =' Y-AXIS  MAGNITUDE  IN  VOLT' 

RETURN 

END 


SUBROUTINE  PREFRE ( TITLE , XS , YS ) 
CHARACTER*64  TITLE, XS.YS 
TITLE=’ SPECTRAL  PLOT' 

XS=' X-AXIS  FREQUENCY  IN  GHZ' 
YS=’ Y-AXIS  MAGNITUDE  IN  VOLT' 
RETURN 
END 


SUBROUTINE  IFFF(N,T,X) 

;*******  ** IrkMekirirkirMrkirkirltMtirkltlrleti  *  *  A-***************  **  *  ***  * 

COMPLEX  W( 1025), X( 1025) 

REAL  T 
INTEGER  MODE 
MODE=l 

CALL  PREFFT( N , MODE , NEXP ,  tf  ) 

CALL  FFT(N,MODE ,T,NEXP,W,X) 

RETURN 

END 


'• ■telr^-kirkiricirk-k-kirk-k'ifk  ******  A*  ft**  ***************  **********  ******  A*  A  A*  Airii  kl 


SUBROUTINE  DECON3 ( K , YSC , Y4 , Y5 , ID , TYPE , DATE , FI , SUBAV , T2 , DATAB ) 

V*  Ar  Jr*  ■ 


COMPUTING  THE  TARGET  RESPONSE  TO  THE  INCIDENT  * 


WAVESHAPE  FOR  THE  COMPUTED  SPHERE.  USING 
RIAD'S  METHOD  OF  OPTIMAL  DECONVOLUTION. 


Xo*Kt(f)=Y5*C0NJ(Y4)*Ysc/(Y4*C0NJ(Y4)+L) 


WHERE 

Xo(f) 


IS  THE  INCIDENT  WAVEFORM 


Ht(  f  ) 


IS  THE  TARGET  TRANSFER  FUNCTION 


Y4 


IS  THE  CALIBRATION  SPHERE 
SUBTRACTED  WAVEFORM  (FREQUENCY) 


Y5 


IS  THE  TARGET  SUBTRACTED  WAVEFORM 
FREQUENCY  DOMAIN 


YSC 


Y1 

Y2 

Y3 

K 

NUMPT 

R 


IS  THE  COMPUTED  SCATTERED  ELECTRIC 
FIELD  FOR  THE  CANONICAL  TARGET 
IS  THE  CALIBRATION  MEASUREMENT 
IS  THE  BACKGROUND  MEASUREMENT 
IS  TTE  TARGET  MEASURED 
IS  THE  NUMBER  OF  POINT 
IS  THE  SAME  AS  K 
IS  THE-  REAL  PART  OF  THE  INVERSE 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


FOURIER  TRANSFORM  OF  Xo*Ht(f)  * 

MAG  IS  THE  MAGNITUDE  OF  Xo*Ht(f)  * 

ID  IS  THE  IDENTIFICATION  STRING  * 

SUBAV  IS  THE  NUMBER  OF  SUB  AVERAGE  * 

DATAB  IS  THE  NUMBER  OF  DATA  BLOCK  * 

FI  IS  THE  DUMMY  STRING  * 

NAME 4  IS  THE  NAME  OF  THE  TARGET  DATA  * 

XS  IS  THE  STRING  FOR  PLOTTING  ROUTINE  * 

FOR  THE  HEADER  * 

YS  IS  THE  SAME  AS  XS  * 

TITLE  IS  THE  TITLE  NAME  FOR  PLOTTING  * 

YN  IS  THE  YES  OR  NO  STRING  * 

L9  IS  THE  SMOOTHING  PARAMETER  * 

TUT  IS  THE  FREQUENCY  DOMAIN  OF  Xo*H(f)  * 

TYPE  IS  THE  WAVEFORM  TYPE  * 

XMI1  IS  THE  MINIMUM  VALUE  * 

XMA1  IS  THE  MAXIMUM  VALUE  * 

DATE  IS  THE  DATE  OF  MEASUREMENT  * 


COMPLEX  Y5( 1025) ,Y4( 1025) ,YSC( 1024) ,TUT( 1024) .SUM 

COMPLEX  L1.L2.L3 

REAL  R( 1024), L9,XMA1,XMI1,MAG( 1024) 

REAL  TIME.T.T2 

INTEGER  K.NUMPT, SUBAV, DATAB 

CHARACTER*16  NAME4 

CHARACTER*64  ID , YS , XS , TITLE .TYPE , DATE , F 1 
CHARACTER*!  YN 


T=T2/FL0AT(K-l)*lE-9 
SUM=(0. 0,0. 0) 

DO  111  1=1, K 

SUM=SUM+Y4( I )*CONJG( Y4( I ) ) 

111  CONTINUE  (  _  _  _ 

WRITE(*!*)  '  ENTER  RELATIVE  SMOOTHING  PARAMETER  :  ' 

READ  (*,*)  L9 
L1=CMPLX(L9 ,0.  0) 

L2=SUM/CMPLX(FL0AT(K+1) ,0.  0)*L1 
DO  1  1=1, K 

TUT( I )=Y5( I )*CONJG( Y4( I ) )*YSC( I ) /( Y4( I )*CONJG( Y4( I ) )+L2 ) 
1  CONTINUE 

TITLE* ' SPECTRAL  PLOT’ 

YS=' DECONVOLUTION' 

XS=' X-AXIS  FREQUENCY  IN  GHZ' 

NUMPT=K 


CALL  MAGNITUD(TUT,NUMPT,MAG) 


XMI 1=0.  0 

XMA1=1.  O/T/2. 0*lE-9 


CALL  P12(TITLE,NUMPT,XMI1,XMA1,MAG,XS,ID,YS,L9) 


CALL  IFFF(K,T,TUT) 


DO  222  1=1, K 
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R(I)=REAL(TUT(I)) 
222  CONTINUE 


XMI 1=0.  0 
XMA1=T2 

XS=' X-AXIS  TIME  IN  NSEC' 
TITLE='TIME  SERIES  PLOT’ 
YS=' DECONVOLUTION’ 


CALL  P12(TITLE,K,XMI1,XMA1,R,XS,ID,YS,L9) 


YN='N’ 

WRITE(*!*)  'WOULD  YOU  LIKE  TO  SAVE  THE  DECONVOLUTION  RESULT  ?: 
READ(*,999)  YN 

IF  (  YN.EQ. 'Y'.  OR.  YN.EQ. 'y')  GO  TO  4376 
GOTO  4388 


4376 

4388 

999 


WRITE(*,*) '  ENTER  FILE  NAME  FOR  THE  TIME  DOMAIN  DATA:  ' 
READ(*,999)  NAME4 

CALL  OUT( NAME4 ,  ID , TYPE , DATE ,  FI , SUBAV ,  T2 , DATAB ,  K ,  R) 

YN='Y' 


WRITE  (*,*)  'REPEAT  WITH  NEW  SMOOTHING  PARAMETER  ?  (Y/N):  ' 
READ(*,999)  YN 

IF  (  YN.  EQ.  'Y' . OR.  YN.  EQ.  'y' )  GO  TO  438 
FORMAT(A) 


RETURN 


END 


C 

C 

C 

C 

C 

c 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


F  =  ARRAY  OF  Y  VARIABLE 
XS=  $  X-SCALE  (STRING) 

YS=  $  Y- SCALE  (STRING) 

ID=  IDENTIFICATION  OF  THIS  GRAPH 

SUBROUTINE  P 1 1( TITLE , NPTS , XMINN , XMAXX , FFF , XS , ID , YS ) 

RM-FORTRAN  Subroutine  to  Plot  a  Solid  Line  with  TITLE 
Based  on  Program  PLOT. FOR.  Dec  1987  by  M. A.  Morgan 

INPUT  DATA  FORMAT 


TITLE  -  64  Space  Title 
NPTS  -  If  Data  Points 
XMIN  -  Real  Min  X  value 
XMAX  -  Real  Max  X  value 
F(N)  -  Input  Data  Array 
JHC  -  Hardcopy  ( l=Yes) 

CHARACTER* 1  YN.DUM 
CHARACTER*64  TITLE, XS, ID, YS 
REAL  X( 1025) ,F( 1025) , XMIN, XMAX, FFF( 1025) , XMINN, XMAXX 
INTEGER*2  N , JROW , JCOL , ISYMBL , I TYPE , CYAN .WHITE , YELLOW 


IHC=0 


WRITE(*,*)  '  WOULD  YOU  LIKE  TO  HAVE  A  HARD  COPY:  ' 

READ(*, 100)  YN 

IF(YN.EQ.  'Y'.OR.YN.EQ.  V)  IHC=1 

WHITE=7 
CYAN=1 1 

DO  1111  1=1 ,NPTS 
F( I)=FFF( I) 

1111  CONTINUE 

XMIN=XMINN 

XMAX=XMAXX 

YELL0W=14 

N=NPTS 

DX=(XMAX-XMIN)/(NPTS-1.  0) 

FMIN=0.  0 
FMAX=0.  0 
DO  22  K=1 ,NPTS 
X(K)=XMIN+(K-1.0)*DX 
IF(F(K). LT.  FMIN)  FMIN=F(K) 

IF(F(K).  GT.  FMAX)  FMAX=F(K) 

22  CONTINUE 

IF(FMIN.GT.  0.  0)  FMIN=0.  0 
IF(FMAX.  LT. 0. 0)  FMAX=0.  0 

C  Computing  Scale  Factors  for  Vertical  Axis 

ABSMIN=ABS(FMIN) 

AB  S  MAX= AB  S ( FMAX ) 

YMAX=AMAX1( ABSMIN , ABSMAX) 

NSCL=INT( L0G10( YMAX) ) 

IF  (YMAX. LT.  1.0)  NSCL=NSCL-1 
YSCL=10.  **NSCL 
FMIN=FMIN/YSCL 
FMAX=FMAX/YSCL 
ABSMIN=ABSMIN/YSCL 
ABSMAX=ABSMAX/YSCL 
DO  33  K=1,NPTS 

33  F(K)=F(K)/YSCL 

YMIN=0.  0 

IF(FMIN. EQ.  0.  0)  GO  TO  37 
35  YMIN=YMIN+0. 5 

IF( ABSMIN. GT. YMIN)  GO  TO  35 
YMIN=YMIN*FMIN/ABSMIN 

37  CONTINUE 

YMAX=0.  0 

IF(FMAX.EQ.  0.  0)  GO  TO  41 
39  YMAX=YMAX+0. 5 

IF( ABSMAX.  GT.  YMAX)  GO  TO  39 
YMAX=YMAX*FMAX/ ABSMAX 

41  CONTINUE 

CALL  QBEEP 

CALL  QSMODE(2) 

WRITE(*,*)  'Press  RET  to  View  Plot.  ..MIT  Again  to  Clear  Screen' 

IFCIHC.NE.  1)  GO  TO  42 

WRITE(*,*)  'A  Hardcopy  Will  Be  Made.  ..Check  that  Printer  is  ON' 

42  CONTINUE 

READ(*, 100)  DUM 

C  Calling  GRAFMATIC  Routines  and  Plotting  Solid  Line  Graph 

ITYPE=1 
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n  c 3000 


ISYMBL=-2 
CALL  QSMODEC 16) 

CALL  QPTXT( 64, TITLE, YELLOW, 16,23) 

CALL  QPTXT( 64 , ID , CYAN , 16 , 19 ) 

CALL  QPTXT( 64, XS, CYAN, 16,20) 

CALL  QPTXT( 64, YS.CYAN, 16,21) 

CALL  QPLOT(110,540,20,240,XMIN,XMAX,YMIN,YMAX,0.  ,0.  ,0,1.  ,1.5) 
CALL  QSETUP(0, CYAN, ISYMBL, CYAN) 

XMAJ0R=XMAX/5.  0 

CALL  QXAXIS(XMIN,XMAX,XMAJOR, 1,1,2) 

YMAJ0R=1.  0 

CALL  QYAXIS(YMIN,YMAX,YMAJOR, 1,1,2) 

JR0W=14+(ABS( YMIN)/(ABS(YMAX)+ABS( YMIN) ) )*220 
JC0L=80-430*XMIN/(XMAX-XMIN) 

CALL  QGTXT(3,'0.  O' .WHITE, JC0L,JR0W,0) 

CALL  QPTXT( 1 , ' s ' .YELLOW, 5 , 16) 

CALL  QPTXT(1, 'c' .YELLOW, 5, 15) 

CALL  QPTXTd,' a' .YELLOW, 5, 14) 

CALL  QPTXTd,'  1'  .YELLOW, 5, 13) 

CALL  QPTXT(1, 'e' .YELLOW, 5, 12) 

CALL  QPTXTd, 'X'  YELLOW, 5, 10) 

CALL  QPTXTC2, 1 10 ' .YELLOW, 4, 8) 

CALL  QPTXT(1, '*' .YELLOW, 5, 7) 

CALL  QPTXT(1,'*' .YELLOW, 5, 6) 


44 

CALL  QCMOV(0,6) 
WRITE(*. 150)  NSCL 
CALL  QTABL( ITYPE,N,X,F) 
READ(*, 100)  DUM 

IF( IHC. NE. 1)  GO  TO  44 
CALL  PRTSC 
CONTINUE 

100 

CALL  QSMODE( 2) 
FORMAT(A) 

150 

FORMAT( 4X, 13) 

RETURN 

END 

C 

C  THIS  SUBROUTINE  PLOTS  TWO  GRAPHS  ON  THE  SAME  AXIS 
C 

C  F  =  Y1  ARRAY  VARIABLE  (  BLUE) 

C  FI  =  Y2  ARRAY  VARIABLE  (GREEN) 

C  AFTER  PLOT  THIS  PROGRAM  BE  AWARE  OF  THE  F1,F 

C  IT  WILL  SENT  BACK  A  SCALE  F,F1 

C 

C  XS  =  $  X-SCALE 

C  YS  ■  $  Y-SCALE 

C  ID  =  IDENTIFICATION 

C*******'ft  ******  ft  ************  ft***  ft  ****'***  ********** 

SUBROUTINE  PL08( TITLE , NPTS , XMlNN.XMAXX , PUM , XS , YS , PUM1 , ID) 


RM-FORTRAN  Subroutine  to  Plot  a  Solid  Line  with  TITLE 
Based  on  Program  PLOT. FOR.  Dec  1987  by  M.  A.  Morgan 

INPUT  DATA  FORMAT 
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TITLE  -  64  Space  Title 
NPTS  -  ft  Data  Points 
XMIN  -  Real  Min  X  value 
XMAX  -  Real  Max  X  value 
F(N)  -  Input  Data  Array 
IHC  -  Hardcopy  (l=Yes) 


CHARACTER* 1  YN.DUM 
CHARACTER*64  TITLE, XS ,ID,YS 

REAL  X( 1025) ,F( 1025) ,F1( 1025) ,PUM( 1025) ,PUM1( 1025) ,XMINN,XMAXX 
INTEGER*2  N,JR0W,JC0L,ISYMBL,ITYPE , CYAN .WHITE , YELLOW , GREEN 

WRITE(*,*)  'WOULD  YOU  LIKE  TO  HAVE  A  HARD  COPY  ?:  ' 

READ(*, 100)  YN 
IHC=0 

IF(YN.EQ. 'Y'.OR.YN.EQ.  V)  IHC-1 
DO  159  1=1, NPTS 
F( I)=PUM( I) 

F1(I)=PUM1(I) 

159  CONTINUE 

XMIN=XMINN 

XMAX=XMAXX 

C******************. 


WHITE=7 

CYAN=11 

YELL0W=14 

GREEN=2 

N=NPTS 

DX=(XMAX-XMIN)/(NPTS-1.  0) 

FMIN=0.  0 
FMAX=0.  0 

FMIN1=0.  0 
FMAX1=0.  0 

DO  22  K=1 ,NPTS 
X(K)=XMIN+(K-1.  0)*DX 
IF(F(K).  LT.  FMIN)  FMIN=F(K) 

IF(F(K). GT. FMAX)  FMAX=F(K) 
IF(F1(K).  LT.  FMIN1)  FMIN1=F1(K) 

IF(F1(K). GT. FMAX1)  FMAX1=F1(K) 

22  CONTINUE 

FMIN=AMIN1(FMIN,FMIN1) 

FMAX=AMAX 1 ( FMAX , FMAX 1 ) 

IF(FMIN. GT.  0.  0)  FMIN=0.  0 
IF(FMAX.LT.  0.0)  FMAX=0.  0 

C  Computing  Scale  Factors  for  Vertical  Axis 

ABSMIN=ABS(FMIN) 

AB  SMAX=AB  S ( FMAX ) 

YMAX=AMAX1( ABSMIN , ABSMAX) 
NSCL=INT( L0G10C YMAX) ) 

IF  (YMAX.  LT.  1.0)  NSCL=NSCL-1 
YSCL=10. **NSCL 
FMIN=FMIN/YSCL 
FMAX=FMAX/YSCL 
ABSMIN=ABSMIN/YSCL 
ABSMAX=ABSMAX/YSCL 
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DO  33  K=1 ,NPTS 

F1(K)=F1(K)/YSCL 

33  F(K)=F(K)/YSCL 

YMIN=0. 0 

IF(FMIN. EQ. 0.  0)  GO  TO  37 
35  YMIN=YMIN+0. 5 

IF(ABSMIN. GT.YMIN)  GO  TO  35 
YMIN=YMIN*FMIN/ABSMIN 

37  CONTINUE 

YMAX=0.  0 

IFCFMAX.EQ.  0.  0)  GO  TO  41 
39  YMAX=YMAX+0. 5 

IF(ABSMAX.  GT.  YMAX)  GO  TO  39 
YMAX=YMAX*FMAX/ABSMAX 

41  CONTINUE 

CALL  QBEEP 

CALL  QSMODE( 2) 

WRITEC*,*)  'Press  RET  to  View  Plot.  ..RET  Again  to  Clear  Screen' 

IF(IHC.NE.  1)  GO  TO  42 

WRITE(*,*)  ’A  Hardcopy  Will  Be  Made.  ..  Check  that  Printer  is  ON1 

42  CONTINUE 

READ(*, 100)  DUM 

C  Calling  GRAFMATIC  Routines  and  Plotting  Solid  Line  Graph 

ITYPE=1 
ISYMBL=-2 
CALL  QSMODEC16) 

CALL  QPTXT( 64, TITLE, YELLOW, 16,23) 

CALL  QPTXT(64 , ID, GREEN ,16,19) 

CALL  QPTXT(64,XS, WHITE, 16,20) 

CALL  QPTXTC 64, YS, CYAN, 16,21) 

CALL  QPLOT( 110,540,20, 240 , XMIN , XMAX , YMIN , YMAX , 0.  ,0. ,0,1.  ,1.5) 
CALL  QSETUP( 0 , CYAN , I SYMBL , CYAN ) 

XMA J 0R=XMAX / 5 .  0 

CALL  QXAXISC XMIN, XMAX, XMAJOR, 1,1,2) 

YMAJ0R=1. 0 

CALL  QYAXIS(YMIN ,YMAX,YMAJOR, 1 ,1,2) 

JR0W=14+( ABS( YMIN)/( ABS( YMAX)+ABS( YMIN) ) )*220 
JCOL=80-430*XMIN/(XMAX'XMIN) 

CALL  QGTXT(3,'0.  0’ .WHITE, JCOL.JROW.O) 

CALL  QPTXTC 1,' s' .YELLOW, 5, 16) 

CALL  QPTXTC 1,'c' , YELLOW, 5 , 15) 

CALL  QPTXTC 1,' a' .YELLOW, 5, 14) 

CALL  QPTXTC 1 , ' 1  * .YELLOW, 5 , 13) 

CALL  QPTXTC 1, 'e' .YELLOW, 5, 12) 

CALL  QPTXTC 1,'X'  YELLOW, 5, 10) 

CALL  QPTXTC 2 , ' 10' , YELLOW, 4, 8) 

CALL  QPTXTC 1,'*' .YELLOW, 5, 7) 

CALL  QPTXTC 1,'*' .YELLOW, 5, 6) 

CALL  QCMOVCO.6) 

WRITEC*, 150)  NSCL 
CALL  QTABLC ITYPE,N,X,F) 

i.  ».  # — I,  I  I..»  I— »— ■  a  a  I  ■«.-!. -■  l  » 

r  »\n  «  »»  jv  n  Jl  h  K  JlTfTf  JTTfTf  TTITn  »  HTfTfTf  TfKTfKH  »  »  A 

C  PLOT  ANOTHER  GRAPH 


GREEN=10 


ISYMBL=-2 


ooononnoooooo  n  n  ci  o  o 


ITYPE=1 

CALL  QSETUP( 0, GREEN, I SYMBL, GREEN) 
CALL  QTABL( ITYPE,N,X,F1) 

RE ADC*, 100)  DUM 
IF(IHC.NE.  1)  GO  TO  44 
CALL  PRTSC 
CONTINUE 
CALL  QSMODE(2) 

FORMAT(A) 

F0RMAT(4X, 13) 

RETURN 

END 


F  =  ARRAY  OF  Y  VARIABLE 
XS=  $  X-SCALE  (STRING) 

YS=  $  Y- SCALE  (STRING) 

ID=  IDENTIFICATION  OF  THIS  GRAPH 

SUBROUTINE  P 1 2 ( TITLE , NPTS , XMINN , XMAXX , FFF , XS , ID , YS , L9 ) 

RM-FORTRAN  Subroutine  to  Plot  a  Solid  Line  with  TITLE 
Based  on  Program  PLOT.  FOR.  Dec  1987  by  M.  A.  Morgan 

INPUT  DATA  FORMAT 

TITLE  “  64  Space  Title 
NPTS  -  #  Data  Points 
XMIN  -  Real  Min  X  value 
XMAX  -  Real  Max  X  value 
F(N)  -  Input  Data  Array 
IHC  -  Hardcopy  (l=Yes) 

CHARACTER*!  YN.DUM 
CHARACTER*64  TITLE, XS, ID, YS.RELAT 
REAL  X( 1025) ,F( 1025) , XMIN, XMAX, FFF( 1025) , XMINN, XMAXX, L9 
INTEGER* 2  N , JROW , JCOL , ISYMBL, ITYPE , CYAN , WHITE, YELLOW 

IHC=0 

WRITE(*,*)  '  WOULD  YOU  LIKE  TO  HAVE  A  HARD  COPY:  ' 

READ(*, 100)  YN 

IF(YN.EQ. 'Y'.  OR.  YN.EQ. 'y' )  IHC-1 

WHITE=7 

CYAN=11 

DO  1111  1=1, NPTS 

F(I)=FFF(T) 

1111  CONTINUE 

.  XMIN=XMINN 

XMAX=XMAXX 
YELLOW=14 
N=NPTS 

DX=(XMAX-XMIN)/(NPTS-1.  0) 

’  FMIN=0. 0 

FMAX=0.  0 
DO  22  K=1 ,NPTS 
X(K)=XMIN+(K-1.0)*DX 
IF(F(K).  LT. FMIN)  FMIN=F(K) 


44 

100 

150 
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IF(F(K).  GT.  FMAX)  FMAX=F(K) 

22  CONTINUE 

IF(FMIN. GT. 0.  0)  FMIN=0.  0 
IF(FMAX.  LT.  0.  0)  FMAX=0.  0 

C  Computing  Scale  Factors  for  Vertical  Axis 

ABSMIN=ABS(FMIN) 

AB  SMAX=ABS ( FMAX ) 

YMAX=AMAX 1 ( ABSMIN , ABSMAX ) 
NSCL=INT( LOG10( YMAX) ) 

IF  (YMAX. LT.  1.0)  NSCL=NSCL-1 
YSCL=10.  **NSCL 
FMIN=FMIN/YSCL 
FMAX=FMAX/YSCL 
ABSMIN=ABSMIN/YSCL 
AB SMAX=AB SMAX/ YSCL 
DO  33  K=1 ,NPTS 


33 

F(K)=F(K)/YSCL 

YMIN=0.  0 

IF(FMIN.  EQ.  0.  0)  GO  TO  37 

35 

YMIN=YMIN+0. 5 

IF(ABSMIN.  GT.  YMIN)  GO  TO  35 
YMIN=YMIN*FMIN/ABSMIN 

37 

CONTINUE 

YMAX=0. 0 

IFCFMAX.  EQ.  0.  0)  GO  TO  41 

39 

YMAX=YMAX+0.  5 

IF(ABSMAX.  GT.  YMAX)  GO  TO  39 
YMAX=YMAX*FMAX/ABSMAX 

41 

WRITE(*,*) 

CONTINUE 

CALL  QBEEP 

CALL  QSMODEC 2) 

'Press  RET  to  View  Plot.  ..RET  Again  to  Clear  Screen 
IF(IHC.NE. 1)  GO  TO  42 

WRITE(*,*) 

'A  Hardcopy  Will  Be  Made. ..  Check  that  Printer  is  ON 

42  CONTINUE 

READ(*, 100)  DUM 

C  Calling  GRAFMATIC  Routines  and  Plotting  Solid  Line  Graph 

ITYPE=1 
ISYMBL=-2 
CALL  QSMODEC 16) 

CALL  QPTXT( 64, TITLE, YELLOW, 16,23) 

CALL  QPTXT(64 , ID, CYAN, 16, 19) 

CALL  QPTXT( 64, XS, CYAN, 16,20) 

CALL  QPTXT( 64, YS, CYAN, 16,21) 

RELAT=' RELATIVE  SMOOTH  =' 

CALL  QPTXT(64,RELAT,CYAN, 16, 18) 

CALL  QQNPUT( 265 , 253 ,L9 , 3) 

CALL  QPLOT(110,540,20,240,XMIN,XMAX,YMIN,YMAX,0.  ,0.  ,0,1.  ,1.5) 
CALL  QSETUP(0, CYAN, ISYMBL, CYAN) 

XMAJ0R=XMAX/5.  0 

CALL  QXAXIS(XMIN,XMAX,XMAJOR, 1, 1 ,2) 

YMAJ0R=1.  0 

CALL  QYAXIS(YMIN,YMAX,YMAJOR, 1,1,2) 

JR0W=14+(ABS( YMIN)/(ABS( YMAX)+ABS( YMIN) ))*220 
JCOL=80-430*XMIN/(XMAX-XMIN) 

CALL  QGTXT(3, ’ 0.0 '.WHITE, JCOL,JR0W,0) 
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CALL  QPTXT(l,'s'  .YELLOW, 5 , 16) 
CALL  QPTXT( 1 , ' c' .YELLOW ,5 , 15) 
CALL  QPTXT( 1 , ' a' .YELLOW ,5,14) 
CALL  QPTXT( 1 , 1 1 ' , YELLOW, 5 , 13) 
CALL  QPTXT( 1, ' e' , YELLOW, 5 , 12) 
CALL  QPTXT(1,'X*  YELLOW, 5, 10) 
CALL  QPTXT(2, ' 10' .YELLOW, 4, 8) 
CALL  QPTXT( 1, .YELLOW, 5, 7) 
CALL  QPTXT( 1, ’*’ .YELLOW, 5, 6) 
CALL  QCM0V(0,6) 

WRITE(* , 150)  NSCL 
CALL  QTABL( ITYPE,N,X,F) 
READ(* , 100)  DUM 
IF(IHC.  NE.  1)  GO  TO  44 


CALL  PRTSC 

44 

CONTINUE 
CALL  QSMODE(2) 

100 

FORMAT(A) 

150 

F0RMAT( 4X , I 3 ) 
RETURN 

END 

SUBROUTINE  PL077(Y1,Y5 ,T2, NUMPT, ID1) 

Y2  =  Y5  =  BACKGROUND  DATA  ARRAY 
Y1  =  SCATTERED  DATA  ARRAY 


* 

* 


22 


INTEGER  NUMPT , DATAB , SUBAV , F 1 , NSHIF , SET 

REAL  Yl( 1024) ,Y2( 1024) ,Y3( 1024) ,Y4( 1024) ,Y5( 1024) 

REAL  XMIN.XMAX 

CHARACTER* 16  NAME1 ,NAME2 ,NAME3,QN 
CHARACTER* 1  YN 

CHARACTER*64  YS.XS, ID, TITLE, DATE 
CHARACTER*64  ID 1, TYPE 1 
CHARACTER*64  ID2.TYPE2 
CHARACTER*64  ID3.TYPE3 
REAL  TIME.T2 

WRITEC*’,*)  ’SHIFTING  SCATTERED  WAVEFORM  &  SUBTRACTING  BACKGROUND' 
WRITEC*,*) 

XMIN=0.  0 
YMAXsT? 

DO  22  K=l, NUMPT 
Y2(K)=Y5(K) 

SET=0 


33  CALL  SHIFTC TITLE, YS.XS, ID) 

QAAA  AA'AAA  AitirA'AitirAirAlrAitirA'AirkrAlrAfrAlrirAirA  A  A  A  A h  it  Jrtrk-kit  it 

CALL  PL08( TITLE , NUMPT , XMIN , XMAX , Y1 , XS , YS , Y2 , ID) 

Q-kitirifirAAirAir/ririflrkirhA-k  AAirirAit-A-Air-IrhirkirfrAirAirkirAAirii  kick  A 

C  SHIFT  THE  CURVE 


CALL  NSHIFT( NSHIF) 

IF(NSHIF.EQ.O)  GO  TO  40 

CALL  GOODCY1.Y2, NSHIF, NUMPT, SET) 

GO  TO  33 
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40  CONTINUE 


SUBTRACT  BACKGROUND  FROM  THE  WAVEFORM  **** 


DO  44  1=1 ,NUMPT 
44  Y3( I)=Y1( I) -Y2( I) 

CALL  SUBTR( TITLE, XS.YS) 

CALL  P11(TITLE,NUMPT,XMIN,XMAX,Y3,XS,ID1,YS) 

WRITEC*,*)  'DO  YOU  WANT  TO  RE-SHIFT  AND  SUBTRACT  ?  (Y/N): 
READ(*,300)  YN 

IF(YN.EQ. 'y’.OR.YN.EQ. 'Y')  GO  TO  33 
DO  77  1=1 ,NUMPT 
77  Y1(I)=Y3(I) 

300  FORMAT(A) 

RETURN 

END 


SUBROUTINE  NSHIFT(NSHIF) 

INTEGER  NSHIF 

WRITEC*  *)  1  ************  *  ***  *****  **********  ****  ******** 1 
WRITEC*!*) 

WRITEC*,*)  ’SCATTERED  WAVEFORM  SHIFTING:  ' 

WRITEC*,*) 

WRITEC*,*)  '+n  shifts  scat  waveform  to  the  RIGHT  n-points' 
WRITE(*,*)  '-n  shifts  scat  waveform  to  the  LEFT  n-points’ 
WRITEC*,*)  ’  0  DISPLAYS  subtracted  waveform  ’ 

WRITEC*,*) 

WRITEC*,*)  ’ENTER  SHIFT  C+/-INTEGER):  ’ 

RE ADC*,*)  NSHIF 

RETURN 

END 


SUBROUTINE  SUBTRC TITLE, XS.YS) 

CHARACTER*64  TITLE, XS, ID, YS 

TITLE=’ SUBTRACTED  WAVEFORM’ 

XS=’ X-AXIS  TIME  IN  NS1 

YS=’ Y-AXIS  MAGNITUDE  IN  VOLT’ 

RETURN 

END 


SUBROUTINE  SHIFTCTITLE,YS,XS,ID) 
CHARACTER*64  TITLE, YS,XS, ID 
TITLE=’TIME  SERIES  PLOT’ 

YS=’ BLUE -TARGET’ 

ID= ’ GREEN - BACKGROUND ’ 

XS=’ SHIFT-TARGET  ONLY* 

RETURN 

END 
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SUBROUTINE  G00D(Y1,Y2, NSHIF, NUMPT, SET) 

* 

PROGRAM  TO  SHIFT  Y2  CURVE  * 

SHIFT  Y1  TO  THE  LEFT  * 

+  SHIFT  Y1  TO  THE  RIGHT  * 

Y1  IS  THE  SHIFTING  CURVE  * 

Y2  IS  THE  BACKGROUND  CURVE  * 

NSHIF  IS  THE  NUMBER  OF  POINT  TO  SHIFT  * 

SET  IS  THE  SUM  OF  THE  POINT  TO  SHIFT  * 

IN  ORDER  TO  ZERO  PADDING  * 

NUMPT  IS  THE  NUMBER  OF  POINT  * 

* 


REAL  Yl( 1024), Y2( 1024), Y3( 1024) 

INTEGER  NSHIF.NUMPT.SET 
5  IF(  NSHIF.  LT.  0  )  THEN 

NSHIF=-1*(NSHIF) 

DO  20  I=NSHIFT , NUMPT-NSHIF-SET 
Y3( I )=Y1( I+NSHIF) 

20  CONTINUE 

DO  10  I=NUMPT -NSHIF -SET, NUMPT 
Y2(I)=0.  0 
Y3( I)=0. 0 

10  CONTINUE 

DO  15  1=1 .NSHIF+SET+l 
Y2( I)=0.  0 
Y3( I)=0.  0 

15  CONTINUE 

ELSE  IF(  NSHIF.  GT.  0  )  THEN 
N=NSHIF+SET+1 
DO  40  I=N, NUMPT 

Y3( I)=Y1( I-NSHIF) 

40  CONTINUE 

DO  30  1=1, N 

Y2(I)=0.  0 
Y3( I)=0.  0 

30  CONTINUE 

DO  35  I=NUMPT-NSHIF-SET, NUMPT 
Y2( I)=0. 0 
Y3(I)=0. 0 

35  CONTINUE 

ELSE 
END  IF 

SET=SET+NSHIF 
DO  50  1=1, NUMPT 
Yl( I)=Y3( I) 

50  CONTINUE 
RETURN 
END 

SUBROUTINE  OUT( NAME , ID , TYPE , DATE , F 1 , SUBAV , T2 , DATAB , NUMPT , Y) 
ID  $  =  CHARACTER  $  TARGET  IDENTIFICATION 
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TYPE  $  =  TYPE  OF  WAVEFORM 

DATE  $  =  DATE  OF  EXPERIMENT  TAKE  PLACE 

FI  =  DUMMY  STRING 

SUBAV  =  SUB  AVERAGE 

NUMPT  =  NUMPT  OF  POINT 

DATAB  =  DATA  BLOCK 

Y  =  Y  ARRAY 

T2  =  TIME  WINDOW 

CHARACTER* 64  ID, TYPE, DATE, FI 

CHARACTER* 16  NAME 

REAL  Y( 1024) ,T2 

INTEGER  SUBAV, NUMPT, DATAB 

0PEN( UN IT= 10 , F I LE=NAME , STATUS® ' UNKNOWN ' ) 

WRITEC10,*)  ID 

WRITE(10,*)  TYPE 

WRITE(10,*)  DATE 

WRITE(10,*)  FI 

WRITE( 10,22)  SUBAV 

WRITE( 10,21)  T2 

WRITE( 10,22)  DATAB 

WRITEC 10,22)  NUMPT 

DO  10  1=1, NUMPT 

WRITE(10,20)  Y(I) 

10  CONTINUE 

20  FORMAT(E12. 6) 

21  FORMATCF5. 2) 

22  FORMAT( 15) 

30  FORMAT(A) 

CLOSEC 10) 

RETURN 

END 
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APPENDIX  B.  THE  ACQUISITION  QUICK  BASIC  PROGRAM 


The  following  computer  program  acquires  the  measurement  data  from  the  HP 
54120T  as  described  in  chapter  5.  The  program  is  written  in  QUICK  BASIC 


1000 

1020 

1040 

1060 

1080 

1240 

1321 

1322 

1324 

1325 

1326 

1327 

1328 

1329 

1330 

1341 

1342 

1343 

1344 

1345 

1346 

1347 
1360 
1380 
1400 
1420 
1440 
1460 
1480 
1500 
1520 
1540 
1560 
1580 
1600 
1620 
1640 
1660 
1680 
1700 
1720 
1740 
1760 
1780 
1800 
io20 
1840 


REM  ***********************  *******  A  A  A  A ■***■*  ***** 


REM  *  MODIFY  4/22/88  * 
REM  *  WRITTEN  BY  LT.  JG  SOONPUEN  SOMAPEE  * 
REM  *  ACQUISITION  * 
REM  *  * 
REM  *  * 


REM  ******************************************* 


D$=DATE? 

DUMMY$=TIME$:  PRINT  "DATE . ";D$:  PRINT  "TIME 

COLOR  7,1 

PRINT  "WELCOME  TO  ACQUISITION  PROGRAM" 

PRINT  "THERE  ARE  3  TYPE  OF  MEASUREMENT" 

PRINT  "1.  BACKGROUND . WAVETYPE" 

PRINT  "2.  CALIBRATE . WAVETYPE" 

PRINT  "3.  TARGET . WAVETYPE" 

PRINT  "  PRESS  THE  FOLLOWING  KEY  " 

PRINT  "  L  (  USED  DATA  LAST  RESPONSE) . 

PRINT  "  D  (  SYSTEM  DEFAULT  ) . 

V$=INKEY$:  IF  V$=""  THEN  1344 
IF  V$="L"  OR  V$  ="1"  THEN  1960 
IF  V$="d"  OR  V$  ="D"  THEN  1360 


GOTO  1344 
REM 

REM  *  DEFAULT 
REM 

OPEN  "DATA. DAT"  FOR  OUTPUT  AS  #1 
Pl=512 
PRINT  #1,P1 


N=1024 


PRINT# 1 ,N 
ND=1 


PRINT#1,ND 

T$="4  inch  Thin  wire" 
PRINT#  1,T$ 

D$=DATE$ 

PRINT# 1,D$ 

T2=l. 563E-10 
PRINT# 1,T2 

W$="Background  Wave  Type" 
PRINT# 1,W$ 

Fl=20 

PRINT# 1, FI 
Sl-16 


PRINT#1,S1 

G1=S1 

^r*T»  A  .  4  II 

i/0£»  y*~~ 

PRINT#1 ,DSE? 


"; DUMMY? 


II 


1860  DUMMY$=TIME$ 

1880  PRINT# 1 , DUMMY  $ 

1900  DELAY 1=38 
1920  PRINT#1,DELAY1 
1940  CLOSE# 1 
1960  REM 

1980  REM  Optionally  include  the  following  declarations 

2000  REM  in  your  program. 

2020  REM  They  provide  appropriate  mnemonics  by  which 

2040  REM  to  reference  commonly  used  values. 

2060  REM  Some  mnemonics  (GETS,  ERR%, 

2080  REM  END%,  ATN%)  are  preceded  by  "B"  in  order  to 

2100  REM  distinguish  them  from 

2120  REM  BASICA  keywords. 

2140  REM 

2160  REM  6PIB  Commands 


2180 

UNL%  =  &H3F 

I 

GPIB  unlisten  command 

2200 

UNT%  =  &H5F 

9 

GPIB  untalk  command 

2220 

GTL%  =  &H1 

1 

GPIB  go  to  local 

2240 

SDOS  =  &H4 

t 

GPIB  selected  device  clear 

2260 

PPC%  =  &H5 

1 

GPIB  parallel  poll  configure 

2280 

BGET%  =  &H8 

V 

GPIB  group  execute  trigger 

2300 

TCT%  =  &H9 

1 

GPIB  take  control 

2320 

LL0%  =  &H11 

9 

GPIB  local  lock  out 

2340 

DCL%  =  &H14 

t 

GPIB  device  clear 

2360 

PPU%  =  &H15 

9 

GPIB  ppoll  unconfigure 

2380 

SPE%  =  &H18 

9 

GPIB  serial  poll  enable 

2400 

SPD%  =  &H19 

1 

GPIB  serial  poll  disable 

2420 

PPE%  =  &H60 

9 

GPIB  parallel  poll  enable 

2440 

PPD%  =  &H70 

9 

GPIB  parallel  poll  disable 

2460 

REM 

2480 

REM  GPIB  status 

bit  vector 

2500 

REM  global  variable  IBSTA%  and  wait  mask 

2520 

BERR%  =  &H8000 

» 

Error  detected 

2540 

TIM0%  =  &H4000 

9 

Timeout 

2560 

BEND%  =  &H2000 

f 

E0I  or  EOS  detected 

2580 

SRQI%  =  &H1000 

9 

SRQ  detected  by  CIC 

2600 

RQS%  =  &HC00 

9 

Device  needs  service 

2620 

CMPL%  =  &H100 

9 

1/0  completed 

2640 

L0K%  =  &H80 

9 

Local  lockout  state 

2660 

REM%  =  &H40 

9 

Remote  state 

2680 

CIC%  =  &H20 

1 

Controller -In -Charge 

2700 

BATN%  =  &H10 

9 

Attention  asserted 

2720 

TACS%  =  &H8 

9 

Talker  active 

2740 

LACS%  =  &H4 

9 

Listener  active 

2760 

DTAS%  =  &H2 

9 

Device  trigger  state 

2780 

DCAS%  =  &H1 

f 

Device  clear  state 

2800 

REM 

2820 

REM  Error  messages  returned  in  global  variable 

2840 

EDVR%  =  0 

1 

DOS  error 

IBERR% 


2860  ECIC%  =  1  '  Function  requires  GPIB-PC  to  be  CIC 

2880  EN0L%  =  2  '  Write  function  detected  no  Listeners 

2900  EADR%  =  3  1  Interface  board  not  addressed  correctly 

2920  EARG%  =  4  '  Invalid  argument  to  function  call 

2940  ESAC%  =  5  '  Function  requires  GPIB-PC  to  be  SAC 

2960  EAB0%  =6  '  1/0  operation  aborted 


2980 

3000 

3020 

3040 

3060 

3080 

3100 

3120 

3140 

3160 

3180 

3200 

3220 

3240 

3260 


3660 

3680 

3700 

3720 

3740 

3760 

3780 

3800 

3820 

3840 

3860 

3880 

3900 

3920 

3940 

3960 

3980 

4000 

4020 

4040 

4060 

4080 


ENEB%  =  7  ' 

E0IP%  =  10  ' 

REM 

ECAP%  *  11  ' 

EFSO%  =12 
EBUS%  =  14  ' 

ESTB%  =  15  ' 

ESRQ%  =  16  ' 

REM 

REM  EOS  mode  bits 
BIN%  =  &H1000  ' 

XE0S%  =  &H800  ’ 
RE0S%  =  &H400  * 


Non-existent  interface  board 
I/O  operation  started  before  previous 
operation  completed 
No  capability  for  operation 
File  system  operation  error 
Command  error  during  device  call 
Serial  poll  status  byte  lost 
SRQ  remains  asserted 


Eight  bit  compare 
Send  E0I  with  EOS  byte 
Terminate  read  on  EOS 


REM 

REM  Timeout  values  and  meanings 


3300 

T10US%  = 

1 

'  Timeout 

of 

10  us  (ideal) 

3320 

T30US%  = 

2 

1  Timeout 

of 

30  us  (ideal) 

3340 

T100US%  = 

3 

'  Timeout 

of 

100  us  (ideal) 

3360 

T300US%  = 

4 

'  Timeout 

of 

300  us  (ideal) 

3380 

T1MS%  = 

5 

1  Timeout 

of 

1  ms  (ideal) 

3400 

T3MS%  = 

6 

'  Timeout 

of 

3  ms  (ideal) 

3420 

T10MSX  = 

7 

1  Timeout 

of 

10  ms  (ideal) 

3440 

T30MS7.  = 

8 

’  Timeout 

of 

30  ms  (ideal) 

3460 

T100MS%  = 

9 

'  Timeout 

of 

100  ms  (ideal) 

3480 

T300MS%  = 

10 

1  Timeout 

of 

300  ms  (ideal) 

3500 

T1S%  = 

11 

'  Timeout 

of 

1  s  (ideal) 

3520 

T3S%  = 

12 

'  Timeout 

of 

3  s  (ideal) 

3540 

T10S%  = 

13 

1  Timeout 

of 

10  s  (ideal) 

3560 

T30S%  = 

14 

1  Timeout 

of 

30  s  (ideal) 

3580 

T100S%  = 

15 

'  Timeout 

of 

100  s  (ideal) 

3600 

T300S%  = 

16 

1  Timeout 

of 

300  s  (ideal) 

3620 

T1000SX  = 

17 

'  Timeout 

of 

1000  s  (maximum) 

3640 

REM 

REM  Miscellaneous 
S%  =  &H8 

LF%  =  &HA  ' 

REM 

REM  Application  program  variables  passed  to 
REM  GPIB  functions 


Parallel  Poll  sense  bit 
Line  feed  character 


REM 

CMD$  =  SPACE$( 10) 

RD$  =  SPACE$( 255) 

WRT$  =  SPACE$( 255) 

BNAME$  =  SPACE$( 7) 

BDNAME$  =  SPACE$( 7) 

FLNAME$  =  SPACE$(50) 

REM 
CLS 

DIM  W( 1500), RRR( 1500) 

DIM  RDSS$( 1024) 

SCREEN  0 
COLOR  7,1,3 

OPEN  "C: DATA. DATA"  FOR  INPUT  AS  #1 
BEEP:  CLS 
REM 


command  buffer 

read  data  buffer 

write  data  buffer 

board  name  buffer 

board  or  device  name  buffer 

file  name  buffer 
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4100 

4120 

4140 

4160 

4180 

4200 

4220 

4240 

4260 

4280 

4300 

4320 

4340 

4360 

4380 

4400 

4420 

4440 

4500 

4520 

4540 

4560 

4580 

4600 

4620 

4640 

4660 

4680 

4700 

4720 

4740 

4760 

4780 

4800 

4820 

4840 

4860 

4880 

4900 

4920 

4940 

4960 

4980 

5000 

5020 

5040 

5060 

5080 

5100 

5120 

5140 

5160 

5180 

5200 

5220 

5240 


REM  *  BEGIN  MAIN  PROGRAM  * 


REM 

LOCATE 

5,8:  PRINT" 

DATA  ACQUISITION  " 

REM 

.  i  - . »  i 

wm 

h**********i 

^** **************************************** 

•Mr 

REM 

* 

T? 

= 

TARGET  IDENTIFICATION 

* 

REM 

* 

w? 

= 

WAVEFORM  TYPE 

k 

REM 

k 

D$ 

= 

DATE 

* 

REM 

* 

N 

2= 

NUMBER  OF  SUB -AVERAGE  /  WAVEFORM  OR  COUNT 

* 

REM 

* 

T2 

= 

DPO  SAMPLING  INTERVAL  (nsec) 

k 

REM 

* 

SI 

= 

DPO  MAXIMUM  DPO  VERTICAL  SCALE(M) 

k 

REM 

k 

ND 

= 

NUMBER  OF  DATA  BLOCKS 

* 

REM 

* 

PI 

= 

NUMBER  OF  POINT  /  WAVEFORM 

* 

REM 

* 

FI 

= 

DPO  TIME  WINDOW 

k 

REM 

* 

DSE$ 

= 

FILE  NAME  OF  THE  DATA  OUTPUT  FILE 

* 

REM 

* 

DEV$ 

= 

NAME  OF  THE  HPDEV  IN  THE  CONFIGURATION  FILE* 

REM 

* 

DSO% 

= 

THE  NUMBER  ASSOCIATED  WITH  THE  DEVICE 

* 

REM 

* 

VEW$ 

= 

STRING?  VARIABLE  FOR  SENDING  TO  HP54120 

* 

REM 

* 

DUMMY $ 

= 

TIME  SET 

* 

REM 

REM 

•frki 

i  ' 

,  n  nA?,  ,i  oTf  AA  M  1 

Hr****  A  A  ****************  A  A  A  A  J.  A  A  A  J.  J.  A  A  A  A  A**** 

** 

INPUT#1,P1 

PRINT  "NUMBER  OF  POINT . PI 

INPUT# 1 ,N 

PRINT  "NUMBER  OF  SUB -AVERAGE ( count ) . ";N 

INPUT# 1,ND 

PRINT  "NUMBER  OF  DATA  BLOCK . ND 

INPUT# 1,T$ 

PRINT  "TARGET  ID . ";T$ 

INPUT#1,D$ 

PRINT  "DATE . D$ 

INPUT# 1 ,  T2 

PRINT  "DPO  SAMPLING  INTERVAL . ";T2 

INPUT# 1,  W$ 

PRINT  "WAVEFORM  TYPE . ";W$ 

INPUT#1 ,F1 

PRINT  "DPO  TIMEWINDOW(  NSEC  ) . ";F1 

INPUT# 1,  SI 

PRINT  "DPO  VERTICAL  SCALE(M) . ";S1 

G1=S1 

INPUT# 1,  DSE$ 

PRINT  "DATA  FILE  NAME  OUTPUT . ";DSE$ 

INPUT#1, DUMMY? 

PRINT  "AUTO  TIME  SET . DUMMY? 

INPUT#1 ,DELAY1 

PRINT  "DELAY  (  NSEC  ) . "; DELAY 1 

CLOSE# 1 

PRINT  " . 

COLOR  7,1:  PRINT  "  PRESS  THE  FOLLOWING  KEY  " 

PRINT  "  Y  (  CHANGE  THE  DATA.  DAT  ) . 

PRINT  "  N  (NO  CHANGE  ) . 

V?=INKEY?: IF  V?=""  THEN  5140 
IF  V?="Y"  OR  V?  ="y"  THEN  11260 
IF  V?="N"  OR  V?  ="n"  THEN  5220 
GOTO  5100 
DEV?="HPDEV" 
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REM  *  OPEN  DEVICE  AND  RETURN  THE  UNIT  * 

REM  *  DESCRIPTOR  ASSOCIATED  WITH  THE  GIVEN  * 
REM  *  NAME  * 

CALL  IBFIND(DEV$,DSO%) 

IF  DS0LU  THEN  PRINT  "ERR9LiN 

REM  »  CLEAR  SPECIFIED  DEVICE  * 

CALL  IBCLR(DSO%) 

GOSUB  6080 
CLS: BEEP 

PRINT  W$+"  DATA  STORAGE  COMPLETE" 

PRINT  "  MEASURE  OTHER  WAVEFORM  TYPE(Y/N)  " 

Q$=INKEY$:  IF  Q$="  "  THEN  GOTO  5560 

IF  Q$="N"  OR  Q$="n"  THEN  GOTO  5640 

IF  Q$="Y"  OR  Q$="y"  THEN  GOTO  4040 

GOTO  5560 

STOP 


REM*  SUBPROGRAM  * 

REM*  * 

REM  *  TEXT  REFERENCE  * 

REM  *  "GURU  II  GPIB  USER'S  RESOURCE"  * 

REM  *  "UTILITY  FOR  IBM  PC"  * 

REM  *  "TEKTRONIX  1986"  * 

REM  *  WRITE  DATA  FROM  STRING  SUBROUTINE  * 

CALL  IBWRT(DS07OJVEW$) 

RETURN 

REM  *  TEXT  REFERENCE  * 

REM  *  "GETTING  STARTED  GUIDE"  * 

REM  *  "HP  54120T  DIGITIZING  "  * 

REM  *  "OSCILLOSCOPE"  * 

REM  *  "HP  COMPANY  1987"  * 

REM  *  SETUP  THE  DISPLAY  * 

VEW$="*RST" 

GOSUB  5840 

REM  *  CRT  FORMAT  1  PROVIDE  ONE  DISPLAY  AREA  * 
REM  *  AND  USE  8  DIVISIONS  FOR  THE  FULL  * 
REM  *  SCALE  RANGE  * 

VEW$=": DISPLAY: GRATICULE  FRAME; FORMAT  l" 

GOSUB  5840 

REM  ***************************  ************** 
REM  *  TURN  TDR  STEP  GENERATOR  ON  * 


WRITE  DATA  FROM  STRING  SUBROUTINE 


TEXT  REFERENCE 
"GETTING  STARTED  GUIDE" 
"HP  54120T  DIGITIZING  " 
"OSCILLOSCOPE" 

"HP  COMPANY  1987" 


SETUP  THE  DISPLAY 


6380 

6400 

6420 

6440 

6460 

6480 

6500 

6520 

6540 

6560 

6580 

6600 

6620 

6640 

6660 

6680 

6700 

6720 

6740 

6760 

6780 

6840 

6860 

6880 

7020 

7040 

7060 

7080 

7100 

7120 

7140 

7160 

7180 

7200 

7220 

7240 

7260 

7280 

7300 

7320 

7340 

7360 

7380 

7400 

7420 

7440 

7460 

7480 

7500 

7520 

7540 

7560 

7580 

7600 

7620 

7640 


VEW$=":  NETWORK:  REFLECTION:  PRESET" 

GOSUB  5840  .  . 

REM  *  PAGE  6-13  TURN  ON  CHANNEL  4 
REM 

VEW$=": VIEW  CHANNEL4" 

GOSUB  5840 
REM 

REM  *  PAGE  6-4  TURN  OFF  CHANNEL!  * 

VEW$=":  BLANK  CHANNELl" 

GOSUB  5840 

*ik,k1e*'jck'fci/rk-k4rk4r/c4(1r/e*-kick**  *  *  irMUr ifrrt  iHHHHHHHr* 

REM  *  SET  UP  CHANNEL  RANGE  AND  OFFSET  * 

REM  *  PAGE  17-6  * 

REM  *  TO  BE  20  NANOSECONDS  * 

REM 

VEW$=":  TIMEBASE:  RANGE  20  NS" 

GOSUB  5840 
REM  PRINT 

REM  PRINT  "SET  THE  VERTICAL  SCALE  RANGE  FROM  8  TO  640  MV 
VERT=S1 

VEW$=": CHANNEL4:  RANGE  "+STR$( VERT)+"m" 

GOSUB  5840 

REM  *  TURN  ON  THE  GRATITUDE  * 

REM 

VEW$=":  DISPLAY:  GRATICULE  GRID" 

GOSUB  5840 

VEW$=":  ACQ:  BAND  LOW" 

GOSUB  5840 
REM 

REM  *  AVERAGING  ON  * 

REM 

VEW$=":  ACQUIRE:  TYPE  AVERAGE" 

GOSUB  5840 
REM 

REM  *  ASKING  FOR  BANDWIDTH  * 

REM 
PRINT 

PRINT  "THE  DEFAULT  BANDWIDTH  IS  12.4GHZ  " 

PRINT  "WOULD  YOU  LIKE  TO  CHANGE  TO  20  GHZ  (Y/N)" 
Q$=INKEY$:  IF  Q$="  "  THEN  GOTO  7380 
IF  Q$="N"  OR  Q$="n"  THEN  GOTO  7500 
IF  Q$="Y"  OR  Q$="y"  THEN  GOTO  7460 
GOTO  7380 

VEW$=": ACQ:  BAND  HIGH" 

GOSUB  5840 
RD$=SPACE$(255) 

VEW$=" ACQUIRE:  COUNT  "+STR$(N) 

GOSUB  5840 

VEW$=": ACQUIRE:  POINTS  "+STR$(P1) 

GOSUB  5840 

REM  *** Irlrlrlrtrkirtrh  ft**  ft  kick  irk***  kit  ********** ft*** 

REM  *  SET  TIME  DELAY  * 

REM  *  * 


I 


\ 


I 


I 

u 


* 


¥ 

| 
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7660 

7680 

7700 

7720 

7740 

7760 

7780 

7800 

7820 

7840 

7860 

7880 

7900 

7920 

7940 

7960 

7980 

8000 

8020 

8040 

8320 

8340 

8360 

8380 

8400 

8420 

8560 

8600 

8620 

8640 

8660 

8680 

8700 

8720 

8740 

8760 

8780 

8800 

8820 

8840 

8860 

8880 

8900 

8920 

8940 

8960 

8980 

9000 

9020 

9040 

9060 

9080 

9100 

9120 

9140 

9160 


REM 
REM 
REM 
REM 

DELAY=DELAY1 
DUM$=STR$( DELAY) 

VEW$=":  TIM:  DEL  r’+DUM$+"NS" 

GOSUB  5840 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 

REM  HP  54120T  SENT  THE  LAST  BYTE  RESPONSE  "LF"  ASII  DECIMAL  10 
VEW$=":  SYSTEM:  HEADER  OFF; : EOI  ON” 

GOSUB  5840 

REM  *********5V***********5'r****************** 

PRINT  "I  AM  GOING  TO  DIGITIZE  CHANNEL(Y/N)" 

Q$=INKEY$:  IF  Q$="  "  THEN  GOTO  8360 
IF  Q$="Y"  OR  Q$="y"  THEN  GOTO  8420 
GOTO  8360 
PRINT  TIME$ 

REM  *  GET  THE  NUMBER  AND  STORE  IN  THE  * 

REM  *  FILE  * 

DSEW$="C:  C" 

FOR  1=1  TO  ND 

IF  I<10  THEN  DSFILE$=DSEW$+RIGHT$(STR$(I),l)+"-WFM" 

IF  I>9  THEN  DSFILE$=DSEW$+RIGHT$(STR$(I),2)+".  VFM" 

REM 
REM 

REM  *  DIGITIZE  CHANNEL  4 
REM 

VEW$=": DIGITIZE  CHANNEL4" 

REM 

GOSUB  5840 

VEW$=": WAVEFORM:  SOURCE  WMEM0RY4;  FORMAT  ASCII" 

GOSUB  5840 

VEW$=":  WAVEFORM:  DATA?" 

GOSUB  5840 
REM 

CALL  IBRDF(DSO%,DSFILE$) 

CLOSE# 1 
NEXT  I 

REM  *  FINISH  DIGITIZING  * 

REM 

PRINT  TIME$ 

DSFILE$="HEADER" 

VEW$=":  WAV:  PRE?" 

GOSUB  5840 
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i 

I 


I 


I 


I 


I 


I 


I 


9180  CALL  IBRDF(DSO%,DSFILE$) 

9200  PRINT  "HEADER  COMPLETE" 

9220  GOSUB  10200 

9240  PRINT  "*********************" 

9260  PRINT  "WOULD  YOU  LIKE  TO  GET  A  HARD  COPY?(Y/N)" 

9280  Q$=INKEY$: IF  Q$="  "  THEN  9280 

9300  IF  Q$="N"  OR  Q$="n"  THEN  GOTO  10080 

9320  IF  Q$="Y"  OR  Q$="y"  THEN  GOTO  9360 

9340  GOTO  9280 

9380  REM  *  PLOT  THE  LAST  WAVEFORM  * 

9400  REM  Jctrklck'k 

9420  VEW$=":  HARD:  SOUR  FACT.WMEM4" 

9440  GOSUB  5840 
9460  PRINT  "PLOT  ?  (Y^Nj)" 


THEN  9480 
THEN  GOTO  9600 


9480  Q$=INKEY$:  IF  Q$=' 

9500  IF  Q$="N"  OR  Q$="n" 

9520  IF  Q$="Y"  OR  Q$="y"  THEN  GOTO  9560 
9540  GOTO  9480 
9560  VEW$=": PLOT?;  " 

9580  GOSUB  5840 
9570  GOTO  9780 

9600  REM  **********************  *  * ***************** 

9620  REM  *  PRINTING  * 

9640  REM 

9660  VEW$=":  HARDCOPY:  PRINTER  DEFAULT" 

9680  GOSUB  5840 

9700  VEW$=": HARDCOPY:  PAGE  AUTOMATIC" 

9720  GOSUB  5840 
9740  VEW$=":  PRINT?;" 

9760  GOSUB  5840 
9780  BOAD$="GPIBO" 

9800  CALL  IBFIND(BOAD$,BOA%) 

9820  CALL  IBSIC(BOA%) 

9840  VEW$="?Q!  " 

9860  CALL  IBCMD  (BOA7.,VEW$) 

9880  V%=1 

9900  CALL  IB GTS  (BOA%,V%) 

9920  PRINT  "■> 

9940  PRINT  "PLEASE  WAIT  UNTIL  THE  PLOTTER  HAD  FINISHED" 
9960  PRINT  "PLEASE  RELEASE  THE  DPO  TO  LOCAL(Y)" 

9980  Q$=INKEY$: IF  Q$="  "  THEN  9280 
10000  IF  Q$="Y"  OR  Q$="y"  THEN  GOTO  10040 
10020  GOTO  9980 
10040  CALL  IBLQC(DSO%) 

10060  CALL  IBLOC( BOA%) 

10080  RETURN 

10100  REM  *  ****************************** ********* 

10120  REM  SUBROUTINE  READ  DATA  TO  STRING  REF.  GURU  D-47 
10140  CALL  IBRD(DSO%,RD$) 

10160  RETURN 

10180  REM  **************  . «*** ****  ***************** 

10200  REM  *  PROGRAM  TO  ARRANGE  FORMAT  FILE  * 

10220  REM  *  SCALE  THE  DATA  * 

10240  REM  *  * 

10260  REM  *  * 
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10300  OPEN  "C:  HEADER"  FOR  INPUT  AS#1 
10320  INPUT#1  ,A%,B%,C%,D/i,E,F,G%,H,  I ,  J,K 
10340  PRINT  "FORMAT  ";A% 

10360  PRINT  "TYPE  B% 

10380  PRINT  "POINT  ";C% 

10400  PRINT  "COUNT  ";D% 

10420  PRINT  "XINCRE  ";E 
10440  PRINT  "XOR  ";F 
10460  PRINT  "XREF  ";G% 

10480  PRINT  "YINC  ";H 
10500  PRINT  "YOR  "j IY 
10520  PRINT  "YREF  ";J 
10540  PRINT  "RANGE  K 
10560  CLOSE# 1 
10580  REM  DSE$="C:C" 

10600  FOR  K=1  TO  ND 

10620  IF  K<10  THEN  DSFILE$=DSEW$+RIGHT$( STR$(K) , 1)+". wfm" 

10640  IF  K>9  THEN  DSFILE$=DSEW$+RIGHT$(STR$(K) ,2)+".  wfm" 

10660  IF  K>99  THEN  DSFILE$=DSEW$+RIGHT$( STR$(K) ,3)+". wfm" 

10680  OPEN  DSFILE$  FOR  INPUT  AS  #1 

10700  REM 

10720  REM 

10740  REM 

10760  REM 

10780  REM 

10800  REM 

10820  FOR  1=1  TO  C?„ 

10840  INPUT  #1,WFM 

10860  W(  I)  •*  ( (WFM-J)*H+IY) 

10880  NEXT  I 

10885  SUMM=0.  0 

10890  FOR  1=1  TO  C% 

10892  SUMM=SUMM+W( I ) 

10896  NEXT  I  _  _ 

10898  REM  *  REMOVE  MEAN  * 

10900  SMEAN=SUMM/C% 

10901  PRINT  "MEAN  ="; SMEAN 

10902  FOR  1=1  TO  C% 

10903  RRR( I)=W( I) -SMEAN 

10904  NEXT  I 
10910  CLOSE #1 

10920  DSFILE$=DSE$+".  WFM" 

10940  OPEN  "A",#1,DSFILE$ 

10960  PRINT# 1 ,T$ 

10980  PRINT#1,W$ 

11000  PRINT# 1,D$ 

11020  PRINT# 1 , DUMMY $ 

11040  PRINT# 1,D% 

11060  PRINT#1, USING  "###.##";  FI 
11080  PRINT#1 ,ND 
11100  PRINT# 1, PI 
11120  FOR  1=1  TO  C% 

11140  PRINT#1, USING"##.#####  ”;RRR(I) 


11160  NEXT  I 
11180  REM 
11200  CLOSE# 1 

11220  NEXT  K 

11221  BEEP: BEEP 
11240  RETURN 
11260  REM  ***** 


11280  REM  *  PI  NUMBER  OF  POINTS  * 

11300  REM  *  N  =  NUMBER  OF  SUBAVERAGES  * 

11320  REM  *  ND  NUMBER  OF  DATA  BLOCKS  * 

11340  REM  *  T?  =  TARGET  IDENTIFICATION  * 

11360  REM  *  D$  DATE  OF  MEASUREMENT  * 

11380  REM  *  T2  —  DPO  SAMPLING  INTERVAL  * 

11400  REM  *  W?  WAVEFORM  TYPE  * 

11420  REM  *  FI  DPO  TIME  WINDOW  * 

11440  REM  *  SI  =  DPO  MAXIMUM  SCALE  * 

11460  REM  *  DSE$  =  OUTPUT  DATA  FILE  NAME  * 

11480  REM  *  DUMMY?  =  DUMMY  ?  * 

11520  TS?=SPACE?( 64) 

11540  DUMMY?=SPACE?(64) 


11560  KEY  OFF:  SCREEN  0,0,0 

11580  OPEN  "DATA.  DAT"  FOR  INPUT  AS  #1 


11600  INPUT#1,P1: INPUT#1,N:  INPUT#1,ND: INPUT# 1,T$:  INPUT#1,D?:  INPUT# 1,T2 
11620  INPUT# 1,W$:  INPUT# 1, FI: INPUT#1,S1:  INPUT# 1 ,DSE?: INPUT# 1, DUMMY? 
11640  INPUT# 1, DELAY 1:  CLOSE# 1 
11660  COLOR  7,1,3:  CLS 
11680  X=4:Y=16 


11700 

11720 

11740 

11760 

11780 

11800 

11820 

11840 

11860 

11880 

11900 

11920 

11940 

11960 

11980 

12000 


PRINT"PROGRAM  DATA  ACQUISTION  " 

LOCATE  X,Y+15:PRINT"PRESS  THE  FOLLOWING  KEY" 

LOCATE  X+2.Y:  PRINT  "  <  F  >  FILE  NAME  FOR  DATA  OUTPUT. . . ! . . 


LOCATE  X+3,Y:  PRINT  " 
LOCATE  X+4,Y:  PRINT  " 
LOCATE  X+5,Y:  PRINT  " 
LOCATE  X+6,Y:  PRINT  " 
LOCATE  X+7,Y:  PRINT  " 
LOCATE  X+8,Y:  PRINT  " 
LOCATE  X+9,Y:  PRINT  " 
LOCATE  X+10,Y:  PRINT  " 
LOCATE  X+11,Y:  PRINT  " 
LOCATE  X+12  ,Y:'PRINT'  " 
LOCATE  X+13,Y:  PRINT  " 
LOCATE  X+14.Y:  PRINT  " 
Q$=INKEY$:  IF  Q?="  "  T 


<  P  >  NUMBER  OF  POINTS . 

<  A  >  NUMBER  OF  SUB  AVERAGES . 

<  B  >  NUMBER  OF  DATA  BLOCKS . 

<  T  >  TARGET  ID.  (?) . ! 

<  D  >  DATE  (?) . 

<  W  >  WAVEFROM  TYPE . ! 

<  K  >  DPO  TIME  WINDOW  . 

"  <  S  >  DPO  SAMPLING  INTERVAL . 

"  <  R  >  RUN  THE  PROGRAM . 

"  <  U  >  DPO  MAXIMUM  VERT  MV . 

"  <  X  >  AUTO  TIME  SET . 

"  <  G  >  DELAY  TIME . 

THEN  12000 


SI 

DUMMY? 

38 


12020  IF  Q?="G"  OR  Q?="g"  THEN  12900 
12040  IF  Q?="F"  OR  Q?="f"  THEN  12300 
12060  IF  Q?="P"  OR  Q?="p"  THEN  12340 
12080  IF  Q?="A"  OR  Q?="a"  THEN  12560 
12100  IF  Q?="B"  OR  Q?="b"  THEN  12600 
12120  IF  Q?="T"  OR  Q?="t"  THEN  12640 
12140  IF  Q?="D"  OR  Q?="d"  THEN  12680 
12160  IF  Q?="W"  OR  Q?="w"  THEN  12700 
12180  IF  Q?="S"  OR  Q?="s"  THEN  12740 
12200  IF  Q?="K"  OR  Q?="k"  THEN  12780 
12220  IF  Q?="R"  OR  Q?="r"  THEN  4040 
12240  IF  Q?="U"  OR  Q?="u"  THEN  12820 
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12260 

12280 

12300 

12320 

12340 

12360 

12380 

12400 

12420 

12440 

12460 

12480 

12500 

12520 

12540 

12560 

12580 

12600 

12620 

12640 

12660 

12680 

12700 

12720 

12740 

12760 

12780 

12800 

12820 

12840 

12860 

12880 

12900 

12920 

12940 

12960 

12980 

13000 

13020 

13040 

13060 


IF  Q$="X"  OR  Q$="x"  THEN  12860 
GOTO  11720 

PRINT" INPUT  FILE  NAME($)  EX. .C:F" 

INPUT  DSE$:  GOTO  12940 

CLS:  LOCATE  4,8:  PRINT"NUMBER  OF  POINT  " 

PRINT  "A).  . .  1024 . POINT.  " 

PRINT  "B) _ 512 . POINT." 

PRINT  "C) _  256 . .POINT." 

PRINT  "D) _ 128 . POINT." 

Q$=INKEY$:  IF  Q$=""  THEN  12440 
IF  Q$="A*  OR  Q$="a"  THEN  Pl=1024 
IF  Q$="B"  OR  Q$="b"  THEN  Pl=512 
IF  Q$="C"  OR  Q$="c"  THEN  Pl=256 
IF  Q$="D"  OR  Q$="d"  THEN  Pl=128 
GOTO  12940 

PRINT" INPUT  NUMBER  OF  SUB -AVERAGE" 

INPUT  N: GOTO  12940 

PRINT  "INPUT  NUMBER  OF  DATA  BLOCK" 

INPUT  ND  : GOTO  12940 
PRINT" INPUT  TARGET  ID($)" 

INPUT  T$:  GOTO  12940 

PRINT"SET  DATE":  D$=DATE$:  GOTO  12940 

PRINT" INPUT  WAVE FROM  TYPE($V 

INPUT  W$:  GOTO  12940 

PRINT" INPUT  DPO  SAMPLING  INTERVAL" 

INPUT  T2:  GOTO  12940 
PRINT"DPO  TIME  WINDOW  " 

INPUT  FI: GOTO  12940 

PRINT" INPUT  DPO  VERTICAL  SCALE" 

INPUT  SI: GOTO  12940 
PRINTOUT 0  TIME" 

DUMMY$=TIME$:  GOTO  12940 
PRINT" INPUT  TIME  DELAY" 

INPUT  DELAY 1:  GOTO  12940 

OPEN  "DATA. DAT"  FOR  OUTPUT  AS#1 

PRINT#1,P1:  PRINT#1,N:  PRINT#1,ND:  PRINT#1,T$:  PRINT#1,D$ 
PRINT# 1 ,T2:  PRINT#1,W$: PRINT#1 .USING  "###. ##";F1: PRINT#1,S1 
PRINT#1 ,DSE$: PRINT#1 .DUMMY? :  PRINT#1 ,DELAY1 
CLOSE# 1  :  CLS 
REM 

GOTO  11700 


APPENDIX  C.  CHARACTERISTICS  OF  LOW  LOSS  CABLE 


Table  21.  CHARACTERISTICS  OF  PORT  RETURN  CABLES 


Impedance 

50  ohms 

Capacitance 

26  pF/ft 

Time  delay 

1.2  ns, 'ft 

Velocity  of  propagation 

85%  of  light  velocity 

Jacket  withstand 

1.0  kV 

Center  conductor 

15  AWG  solid 

Minimum  bend  radius 

1  inch 

Dielectric  constant 

1.4 

Dielectric  withstand 

1.0  kV 

Table  22.  INSERTION  LOSS 


0-4  GHz 

less  than  0.36  dB 

4-8  GHz 

less  than  0.5  dB 

8-12  GHz 

less  than  0.6  dB 

12-16  GHz 

less  than  0.72  dB 

16-18  GHz 

less  than  0.8  dB 

18-26.5  GHz 

less  than  1.0  dB 
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